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LIMITING DISTRIBUTION OF ELLIPTIC HOMOGENIZATION ERROR 
WITH PERIODIC DIFFUSION AND RANDOM POTENTIAL 

WENJIA JING 


Abstract. We study the limiting probability distribution of the homogenization error for second 
order elliptic equations in divergence form with highly oscillatory periodic conductivity coefficients 
and highly oscillatory stochastic potential. The effective conductivity coefficients are the same as 
those of the standard periodic homogenization, and the effective potential is given by the mean. 
We show that in the short range correlation setting, the limiting distribution of the random part 
of the homogenization error, as random elements in proper Hilbert spaces, is Gaussian and can 
be characterized by the homogenized Green’s function, the homogenized solution and the statistics 
of the random potential. Similar results hold for random potentials that are functions of long 
range correlated Gaussian random fields. These results generalize previous ones in the setting with 
slowly varying diffusion coefficients, and the current setting with fast oscillations in the differential 
operator requires new methods to prove compactness of the probability distributions of the random 
fluctuation. 

Key words: periodic and stochastic homogenization, random field, probability measures on Hilbert 
space, weak convergence of probability distributions 


1. Introduction 


In this article we study the limiting distribution, in certain Hilbert spaces, of the homogenization 
error for second order elliptic equations in divergence form with highly oscillatory periodic diffusion 
coefficients and highly oscillatory random potential. 

More precisely, we consider the following Dirichlet problem on an open bounded subset D C M n , 
with homogeneous boundary condition and a source term / € L 2 (D), 

(Mj (7) ^{x,u)\ + ? (f’ w ) “ £ ( x ’ w ) = /( x )> 


dxi 


x £ D, 
x € dD. 


( 1 . 1 ) 


u £ {x) = 0, 

The conductivity coefficients (ajj(-)) and the potential are highly oscillatory in space, and 

0 < e <C 1 indicates the small scale on which these coefficients oscillate. We assume that the 
conductivity coefficients are deterministic and periodic, and the potential is a stationary random 
field on some probability space P). More precise assumptions are given in section [2j It is 

well known that, under mild assumptions like stationary ergodicity of q(x, lo), the equation above 
homogenizes, i.e. u £ converges, almost surely in H, weakly in H l (D) and strongly in L 2 (D ) to the 
solution of the deterministic homogenized problem 

d 2 u 

'«a^fc“ (l)+5 " (l) = /(l) ' x€D - (1,2) 


&o 


u(x) = 0, 


x € dD. 
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Here, the effective conductivity coefficients {a t j) are constants defined by 



(1.3) 


where T d = [0, \} d denotes the unit cell and the correctors x k i k = 1, ■ ■ • ,d, are given by the unique 
solution of the corrector equation 



with the normalization condition J Td x k dy = 0; e*, above is the kth standard unit basis vector of 
M d . We note that this formula for (a^) is exactly the classic periodic homogenization formula for 
effective conductivity. The effective potential q in (11.21) is given by the constant 


q = Eq(0,u), 


(1.5) 


where E denotes the mathematical mean with respect to P. 

In this paper we study the law (probability distribution) of the homogenization error u £ — u, 
viewed as random elements in certain Hilbert spaces. We split this error as two parts: E u £ — u 
and u £ — E u £ . In view of the deterministic oscillations in the diffusion coefficients, we expect 
that the periodic homogenization error, in the replacement of (<W( A )) to (dij), makes significant 
contributions to the deterministic error Eit £ — u. Indeed, we show later that this error is essentially 
of order O(e), the same as periodic homogenization. On the other hand, the effect of the random 
potential q(^,uj) becomes visible in the random fluctuation u £ — E u £ , in which the (large) mean is 
removed. We are interested in characterizing the size and the law of this random fluctuation, and 
the answers depend on finer information of the random potential q, such as the decay rate of the 
correlations in q and higher order moments of q ; see section [2] for notions and definitions. 

We find that, when q(x, oj) has short range correlations, the random fluctuation u £ —Eu £ scales like 
£ 2 a2 in L 1 (H,L 2 (Z1)) norm, and scales like £2 when integrated against a test function. Moreover, 
the law of the scaled random fluctuation e~^{u £ — Eii £ ) in L 2 (D) for d = 2,3 and in H~ 1 (D) for 
d = 4,5, converges to Gaussian distributions as follows (see Theorem 12.41 for details): 



Here, W(y ) is the standard multiparameter Wiener process and hence the law of the right hand 
above defines a Gaussian probability measure on L 2 (D) or i7 _1 (D). This Gaussian distribution is 
determined by: G(x,y) which is the Green’s function associated to the homogenized problem (11.21) . 
u{y) which is the homogenized solution, and a which is some statistical parameter of the random 
potential q(x,u). 


We also consider the case when q(x,u) has long range correlations and is constructed as a 
function of Gaussian random field. Then the random fluctuation scales like £" 2 A2 in L 1 (H,L 2 (T>)) 
norm, and scales like £"2 when integrated against a test function. Moreover, the law of the scaled 
random fluctuation s~ "2 {u £ — Eu £ ) in L 2 (D) for d = 2,3 and in H^ 1 (D) for d = 4, 5, converges to a 
Gaussian distribution that can be written as a stochastic integral as above, but with dW replaced 
by W a dy where W° is a centered Gaussian random held with correlation function \x — y | - “, and 
a replaced by some other statistical parameter; see Theorem 16.21 for details. 
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The study of limiting distribution of the homogenization error goes back to m, where the 
Laplace operator with a random potential formed by Poisson bumps was considered. General 
random potential with short range correlations was considered recently in [4], and in [3 |S] for 
other non-oscillatory differential operators with random potential. Long range correlated random 
potential was considered in [5|. When oscillatory differential operators are considered, limiting 
distribution of homogenization error was obtained in m for short range correlated elliptic coeffi¬ 
cients, and in [6] for long range correlated case, all in the one dimensional setting. The main results 
of this paper show that the general framework developed in [4] [8, [5], in order to characterize the 
random fluctuation caused by the random potential, applies even when there are oscillations in 
the differential operators, as long as these oscillations are not statistically related to those of the 
random potential. 

Our approach is as follows: we introduce an auxiliary problem with periodic diffusion coefficients 
and homogenized potential; let v £ be the solution. Then the deterministic homogenization error 
Eu £ — u is essentially characterized by v £ — u, which amounts to classical periodic homogenization 
theory. The random fluctuation u £ — E u £ is then the same as ( u £ — v £ ) — E(n £ — v £ ), which can 
be represented as a truncated Neumann series. The first term X £ in this series contributes to the 
limiting distribution. By Prohorov’s theorem, we need to show that the probability measures of 
{X £ } is tight in the proper Hilbert space, and that their characteristic functions converges. The 
latter is essentially the convergence in distribution of the integration of X £ against test functions, 
in view of the uniform in e estimates of the Green’s functions associated to the oscillatory diffusion, 
this step is the same as the earlier setting with non-oscillatory diffusion. The role of oscillations in 
the diffusion, however, becomes prominent in the step of proving tightness of the measures of {X £ }. 
The simple and natural method used in |5j fails completely; see section [3 for details. New ideas are 
needed: we obtain tightness of the the measures of {X £ } in L 2 (D ) by controlling the mean square 
of H s norm of {X £ }, for some 0 < s < similarly, we get tightness in H^ 1 (D) by controlling 
the mean square of H~ s norm with | < s < 1 . The constraints on the spatial dimension d arises 
naturally in the proof of such controls. 

Our analysis relies on uniform estimates of the Green’s function associated to the periodic homog¬ 
enization problem; we refer to mm for the classical results, and to [20, [21] for recent development 
in this direction. We refer to Eauu for recent results on uniform estimates of Green’s function 
for equations with highly oscillatory random diffusion coefficients in spatial dimension higher than 
one. We remark also that in this setting, limiting distribution of the corrector function and that 
of the full random fluctuation u £ — E u £ , in negative Holder space, were obtained in [26j and m 
respectively, in the discrete setting; see also E3- Such results are apparently more challenging to 
obtain, and the proofs are based on calculus in the (infinite dimensional) probability space. 

The rest of this paper is organized as follows: in section [2] we make precise the main assumptions 
on the parameters of the homogenization problem, in particular on the properties of the random 
potential, and state the main results in the short range correlation setting. Homogenization of (11.ip 
and some useful results on periodic homogenization theory are recalled in section [3] Section [4] and 
[5] are devoted to the proof of the main results, where we characterize how the random fluctuation 
scales in energy norm and in the weak topology, and determine the limiting distribution of the 
scaled fluctuation. We present new methods to prove the tightness of the probability measures 
of the random fluctuations. In section [6] we state and prove the corresponding results in the long 
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range correlation setting. We make some comments and further discussions in section [7] and prove 
some technical results, such as tightness criterion for probability measures, in the appendix. 


2. Assumptions, Preliminaries and Main Results 

2.1. Assumptions on the coefficients. Throughout this paper, we assume that the domain D 
in (11.11) is an open bounded set of R d with C 1,1 boundary. The coefficients cuj(|) and are 

the scaled version of a,ij(x ) and q(x,u). We make the following main assumptions on a^j and q. 

Periodic diffusion coefficients. For the functions (a^), we assume: 

(Al) (Periodicity) The function A := (ctij) : —»• R dxd is periodic. That is, for all x G R d , 

k € 7L d and i,j = 1 , 2 , • • • , d, we have 

aij{x + k) = aij(x). (2.1) 

(A2) (Uniform ellipticity) For all y G T rf , the matrix A(y) = ( a,ij(y )) is uniformly elliptic in the 
sense that, for all £ G R d , one has 

d 

A|£| 2 < Z T A(y)Z = J2 Zi"ij(uKi < Aiei 2 (2-2) 

*j'=i 

(A3) (Smoothness) For some 7 , M with 7 € (0,1] and M > 0, one has 

||^4|lcT'(T d ) ^ At. (2.3) 

We henceforth refer to the above assumptions together as (A). 

Random potential. For the random held q(x,u) on the probability space (U, J 7 , P), we assume: 

(P) (Stationarity and Ergodicity) There exists an ergodic group of P-preserving transformations 
(r x ) 3 , gRt i on U, where ergodicity means that E € T and 

t x E = E, for all x G R. d 

imply that P(E) € {0,1}. The random potential q(y,u) is given by q(T y uj) where q : U —>• R 
is a random variable satisfying 

0 < q(uj) < M, for all w £ fl. (2.4) 

Further assumptions on q. The above assumptions are sufficient for proving homogenization 
result. However, to estimate the size of the homogenization error and to characterize the limiting 
distribution of the random fluctuation, more assumptions on the random held q(-,u) are necessary. 
To simplify notations, we write in the sequel 

q(x,u>) = q + u(x,ui) 

where q is the mean of q and v is the huctuation. Note that q is a deterministic constant and v is 
a mean zero stationary ergodic random held. The auto-correlation function R(x) of q (and hence 
u) is defined as 

R(x) = E [v(x + y, u)v(y, w)] , a 2 := / R(x)dx. (2.5) 

J R d 
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By Bochner’s theorem, R(x) is a positive definite function and a 2 > 0. We assume that a > 0. 
When R is integrable on M d , i.e. a 2 < oo, we say that q has short range correlations; we say q has 
long range correlations if otherwise. 

Short range correlated random fields. In this case, we make an assumption on the rate of 
decay of the correlation function. We denote by C the set of compact sets in W l , and for two sets 
K ], K ‘2 in C, the distance d(K \, K 2 ) is defined to be 

d(Ki,K 2 )= min \x — y\. 
xeK 1 ,yeK 2 

Given any compact set K C C, we denote by J~k the cr-algebra generated by the random variables 
{q(x) : x E A"}. We define the “maximal correlation coefficient” g of q as follows: for each r > 0, 
g(r) is the smallest value such that the bound 

E(¥>i(g)¥> 2 (g)) < g{r)yj'E.((p{{q)) E (tp%(q)) (2.6) 

holds for any two compact sets K 1 , K 2 E C such that d(K\, K 2 ) > r and for any two random 
variables of the form tpi(q), i = 1,2, such that <fi(q) is -measurable and E ipi(q) = 0. We assume 
that 

(S) The maximal correlation function satisfies g 2 E A 1 (M_ t _, r d ~ 1 dr), that is 

g 2 ( r)r d ~ 1 dr < 00 . 

Assumptions on the mixing coefficient g of random media have been used in [3J I8j [IB] ; we refer 
to these papers for explicit examples of random fields satisfying the assumptions. We note that the 
auto-correlation function R(x) can be bounded by g. For any x € M . d , 

|A(x)| = |E(g(s) - Eq)(q(0) - Eg)| < e (|s|)Var( 9 ). 

By (12.41) . q and hence its variance is bounded. In view of (S) and the fact that one can assume 
g E [0,1] (hence g < y/g), we hnd that R is integrable. Therefore, (S) implies that q(x,u) has 
short range correlations. In fact, (S) is a much stronger assumption, and not necessary for the 
main results of this paper to hold. In section [7] we will provide alternative and less restrictive 
assumptions that is sufficient. However, using the assumption (S) and Lemma 14.31 below, we can 
simplify significantly certain fourth order moment estimates of the random potential u(x,u ;); such 
estimates appear often in the study of limiting distribution of the homogenization error. 

Long range correlated potentials constructed from Gaussian fields. There is no general 
tools like central limit theorems when the random held has long range correlations. We restrict to 
a specific family of random fields, constructed from long range correlated Gaussian random fields. 
Let q(x,u) =q + v^x,uj) with q being a nonnegative constant; we assume 

(LI) u(x, oj) = &(g(x)); g(x, u) is a centered stationary Gaussian random held with unit variance. 
Furthermore, the correlation function of g(x,u) has heavy tail. That is, for some positive 
constant n g and some real number a E (0, d), 

R g (x) := E {g(y, u})g(y + x, cj)} ~ K g \x\~ a as |x| —>• 00 . 
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(2.7) 



( 2 . 8 ) 


(L2) The function $ : R —>• M satisfies —q < $ < M — q, and has Hermite rank one, i.e. 

I &(s)e~^~ds = 0 , V\ := J s&(s)e~^~ds / 0 . 

Jr Jr 

(L3) The Fourier transform $ of the function <L satisfies 

/ l^>(Ol(l + lei 3 ) < oo. (2.9) 

JR 

We henceforth refer assumption (L) to the above conditions all together. 

The assumption (L2) makes v(x,lu) = &(g(x,u)) mean zero, and the bounds on <h ensures that 
0 < q(x,oj) < M, which is (12.41) . From the above construction, we check that u(x,uj) is stationary 
ergodic and has long-range correlation function that decays like /e|a;| -a , n = V 2 K g \ see Lemma lA.31 
for the details. Assumption (L3) allows one to derive a (non-asymptotic) estimate, see Lemma I A. 5 1 
in the appendix, for the fourth-order moments of v{x,oo,uj). 

Notations. Throughout the paper, by universal parameters we refer to A, A, 7 and M in the 
assumptions (A), the autocorrelation function R, a 2 , and the mixing coefficients g in the short 
range correlation setting, a, R g , and n in the long range correlation setting, the domain D 

and its boundary dD, and the dimension d. If a constant C depends only on these parameters, 
we say either C depends on universal parameters or C is a universal constant. For the random 
potential v(x,ui) and the functions g(x), R(x), etc., which are related to v. we use is £ ,g £ ,R £ , etc., 
to denote the scaled version. For instance, n £ (x,oj) is short-hand notation for i/(|). We use the 
notation H S (K ), s > 0, for the Sobolev or the fractional Sobolev space W S,2 (K) on some domain 
K C M . d ; when K is bounded, we use Hq(K) for the subspace that consists of functions having trace 
zero at dK ; note that Lfg(M d ) = H s (R d ). We denote by H~ S (K), s > 0, the dual space (Hq(K))'. 
For any Hilbert space R, (•, •)% denotes the inner product in R; when R = L 2 (D), we very often 
omit the subscript and write (•, ■) instead. We use (/, g) whenever the formal integral f D fg makes 
sense. We typically use 1a for the indication function of a set A C M d , or if A is a statement the 
indication function of A being true. Finally, for two real numbers a and b, a A b is a short-hand 
notation for rninja, b}, and a V b means max{a, b}. 

2.2. Probability distribution on functional spaces. We view the random fluctuation u £ — Ert £ 
in the homogenization error as random elements in certain functional spaces, and aim to find the 
limit of its law in that space. It turns out that the choice of functional spaces depends on the 
spatial dimension d. 

When d = 1, one can choose the space C(D ) of continuous functions. In fact, convergence 
in distribution in C(D ) was proved in |4], for random diffusion coefficient a(x,uj) with random 
potential q(x, u), both having short range correlations. In this paper, we prove that for d = 2,3, 
the space can be chosen as L 2 (D ) and for d = 4,5, the space can be chosen as H^ 1 (D). Note that 
both choices are Hilbert spaces. We recall some facts concerning weak convergence of probability 
measures on Hilbert spaces. We refer to the books of Billingsley [9] and Parthasarathy [ 29] for 
more details. 

Probability distributions on a Hilbert space. Let R be a separable Hilbert space, and let X(u) be 
a R valued random element on the probability space (fl, J~, P). Then X determines a probability 
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measure P x on where B(LL) denotes the Borel a- algebra generated by open sets in TL, 

by 

P X {S) = F(X E S), for any S G B(H). (2.10) 

We say a family {^f e } £ g(o,i) of random elements in 7i converge in probability distribution (or in 
law), as e —>• 0, to another random element X on if the probability measures P x ‘ converges 
weakly to P x , i.e. for any real bounded continuous functional / : LL —> R, 

[ f(g) dP x \g) —► [ f(g) dP x (g). 

Jn jh 

In particular, any probability measure P on a separable Hilbert space Li is determined by its 
characteristic function cj) P : LL —>• C, 

<f> P (h) = [ e^* dP(g), (2.11) 

JH 

Moreover, the following result holds: 

Theorem 2.1 ( |29l Chapter VI, Lemma 2.1]). Let {V e } eg (o,i) an d X be random elements in Li, 
possibly defined on different probability spaces. X £ converges to X in law in Li, as e ^ 0, if the 
family of probability measures {P Ae } eg ( 0 i i) tight and for any h E Li, 

lim cj) P (h) = (j) P ( h ). ( 2 - 12 ) 

£->0 

Remark 2.2. Let H = L 2 (D), which is a separable Hilbert space, and let X be a random element 
in L 2 (D ) defined on the probability space P). The characteristic function of P x can be 

calculated as follows: for any h € L 2 (D), 

fi pX {h) = [ e iz dP x ({(h, g) > z}) = [ e iz d¥({(h, X(u)) > z}) =Ee i( ^. (2.13) 

J r Jr 

Therefore, to prove that X s converges in distribution to X as L 2 paths, it suffices to show that 
{P A£ } is tight and that for any h G L 2 (D ), 

{h, X £ ) -> (h, X), (2.14) 

that is, the random variables (h,X £ ) converges in distribution to the random variable (. h,X ). 

In Theorem lA.il in the appendix, we provide a tightness criterion for {P A£ } on L 2 (D), with the 
assumption that {V E (-,u;)} are in Hq(D) for certain s > 0. The criterion is sufficient but by no 
means necessary. Nevertheless, it is very handy for our analysis since the random fields X £ that 
we are dealing with come from solutions of ( 11 . 11 ) . and hence are naturally in Hq(D). 

2.3. Main results. We now state the main results of the paper under the assumption that q(x,u ) 
has short range correlations. Analogous results for long range correlation setting will be presented 
in section [ 6 l 

The first main theorem concerns how does the homogenization error scales. 

Theorem 2.3. Let D C be an open bounded C 1,1 domain, u £ and u be the solutions to (11.11) 
and m respectively. Suppose that (A)(P) and (S) hold, f € L 2 {D) and 2 < d < 7. Then, there 
exists positive constant C, depending only on the universal parameters, such that 

^\W £ ~u \\ L 2 < Ce\\f\\ L 2 . 
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(2.15) 





Moreover, 


E \\u £ — Ku £ 


C'^II/IIlL if 4, 

Ce 2 \ loge|5 ||/|| L 2 , if d = 4. 


(2.16) 


II 2 < 


Moreover, for any ip € L 2 (D), 


E |(u e — Eu e , </?) L 2| < Ce 2 


L 2 11/ IIl 2 


(2.17) 


This theorem provides L 1 (f2, L 2 {D)) estimates of u £ — u and its random part, and its proof is 
detailed in section [4| We note that the size of the full homogenization error is much larger than 
that of its random part. This is because the oscillations in the diffusion coefficients cause some 
deterministic fluctuation in the solution of size 0{e), as in standard periodic homogenization. The 
additional random fluctuation caused by the short range correlated random potential scales like 
£ (dA4)/2 j n ener gy norm, and scales like e d / 2 in the weak topology. These results agree with the case 
of non-oscillatory diffusion coefficients; see @1 IS]. The next result exhibits the limiting law of the 
rescaled random fluctuation e~^{u £ — E u e ). 

Theorem 2.4. Suppose that the assumptions in Theorem \2.3\ hold. Let cr be defined as in (E3|) 
and G(x,y ) be the Green’s function of fll,2|> . Let W(y ) denote the standard d-parameter Wiener 
process. Then 

(i) For d = 2, 3, as e —\ 0, 

'MML „ [ G(x,y)u(y)dW(y), in L\D). (2.18) 

V £ d Jd 

(ii) For d = 4,5, as e —> 0, the above holds as convergence in law in H^ 1 (D). 


The proof of item (i) above can be found towards the end of section 15.11 and that of item (ii) is 
at the end of section IT 21 

Remark 2.5. The integral on the right hand side of (12.181) is understood, for each fixed x, as 
Wiener integral in y with respect to the multiparameter Wiener process W{y). Let X denote the 
result. For d = 2,3, because the Green’s function G(x, y) is square integrable, X is a random 
element in L 2 {D). For d = 4, 5, X is understood through the Fourier transform of its distribution: 
given h* € (j> pX (h*) is defined to be Eexp (ia f D (G(-,y), h* (-))u(y)dW(y)), where E is 

the expectation with respect to the law of W. 


3. Homogenization and Periodic Error Estimates 

The following homogenization result for (11.11) . without the random potential q £ (x, u), is well 
known. The effect of the presence of q £ turns out to be minor for homogenization; nevertheless, we 
include a proof here for the sake of completeness. 

Theorem 3.1. Assume (A1)(A2) and (P). Then there exists Hi € T such that P(Hi) = 1, and 
for all uj £ Hi, the solution u £ of CD converges to the solution u of m weakly in H l {D) and 
strongly in L 2 (D), for any f £ H~ l {D). 
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+ q, 


Let C £ denote the differential operator 

_ _d_ 

dxi 


(3.1) 






and let C £,u be the differential operator C £ + We remark that C £ has highly oscillatory 

but deterministic coefficients while C £,UJ has, in addition, a highly oscillatory and random potential. 
Let Q £,lJ and Q e be the solution operator of the Dirichlet boundary problems associated to C £,u> and 
C £ . Owing to the conditions (12.21) and (12.41) . Q e,w is well defined for any Moreover, we have 

the standard estimate, for any uj G and e > 0, 

\\G £,lxl f\\H 1 (D) < C\\f\\ H -i( D ) (3.2) 

with some constant C that depends on on the universal parameters, and neither on ui nor e. By 
the same token, Q e is well defined and shares the same estimate above. 

Proof of Theorem \3.1\ Step 1. For each uo G 0, the solution u £ of (11.11) is given by C £ ' u /, which 
satisfies the standard estimates 

ll^llz/RD) + ||^4 e V'U e ||^2(£)) + \\q‘-(x,u})u £ \\ L 2 ( D ' ) < C, 

where C depends the universal parameters and /, and is uniform in e and co. As a result, due to 
the compact embeddings H 1 (D) ^ L 2 (D) through a subsequence £j{u) —>• 0 which 

by an abuse of notation is still denoted by e, we have 

Vix E (-, uj) V«(-,«), a(-) Vu £ (-,u) V£(-, W ), 

£“►0 \£/ 6 ~>0 /q o \ 

L 2 / • \ jj-i ' ’ ' 

u £ (-,u) —q [-,(*}) u £ (-,uj) — 

£—^0 \£ J £—>0 

for some function v(-,oj) G H l (D) and some vector valued function £(■,<*;) G [L 2 (D)] d . 

Step 2. Recall that {x k } k =\ are the correctors defined in (|1.4|) . and we can extend them period¬ 
ically to functions defined on M d . Since A(y)(ek + V\ ,fc (y)) is periodic, we have that 

^(e) ( 6fc+ () (f)) J Td A (y">( e k(y) + ^X k (y))dy = Ae k . (3.4) 

For the same reason and the fact f Jd Vx k dy = 0, we have 

e-k + ( v X fc ) (^) J e k + V x k (y)dy = e k . (3.5) 

Now £x an arbitrary function <p G Cq°(D). For each fixed ui G O, let e{uj) —>• 0 be the subsequence 
in Step 1. Consider the integral 

J A (^j X7u £ (x, uj) ■ V jx*; + ex k (-y) | <p(x)dx. 

On one hand, in view of the third item in (13.31) . (13.51) . and the facts that div(A e Vn £ ) = —/ + q £ u £ 
converges in (to -f+p(-,oj) where p is defined in (13.31) 1 and that e k + (Vy fc )(a:/e) is curl-free, 
by the div-curl lemma HI Lemma 1.1], the above integral satisfies 

lim J A X7u £ (x,u>) ■ V ^x k + (!)} c p(x)dx = J f(x,oj) ■ e k ip{x)dx. 

On the other hand, in view of the first item in (13.31) . (13.41) . and the facts that div(yl £ (efc-|-V / \ fc (x/e))) 
converges in H (they all equal to zero) and that \7u £ is curl-free, by the div-curl lemma, we have 

lim Id A^Vu E (x,uj)-v{x k + £x k (!)} <p(x)dx = J X7vA ■ e k <p(x)dx. 

The two limits above must be equal, and it follows that £(•,<*;) = AVv(-,uj) in distribution. 
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Step 3. Recall that the stationary random potential q(x,ui) can be written as q(r x uj) where q 
is an essentially bounded random variable on fi. By Birkhoff ergodic theorem m Theorem 7.2], 
there exists fii E T with P(fii) = 1, and for each wEfii, the whole sequence 


q (p w) =q I ' - 1 ° - c(R - q = Eg(0, u), (3.6) 

for any a £ (l,oo). From the weak formulation of u £ , for any w E 11 1 and for any cp £ Cfi°(D), we 
have 

j A ^ \7u £ (x,u)'V(p(x)dx + j q( — ,u?j u £ (x)<p(x)dx = j f(x)ip(x)dx. 

Pass to the limit along the subsequence e[uj) found in Step 1, we have 


/ A'Vv-'V<p + / qv(x)ip(x)dx = / f(x)<p(x)dx — lim / q ( —, cj) (u £ — v)ip(x)dx. 

Jd Jd Jd Jd ' 

The first term on the left follows from the div-curl lemma and the fact that £ = RVu; the second 
term on the left is due to (13.61) . Finally, the last term on the right hand side is zero since q is 
uniformly bounded and u £ — v converges to zero strongly in L 2 (D). Consequently, the above limit 
shows that v solves the homogenized equation (11.21) . By uniqueness of the homogenized problem, 
we must have that v = u and v is deterministic. 


Finally, for each u £ fij, by the weak compactness in H 1 (D) and the uniqueness of the possible 
limit, the whole sequence u £ converges to u. This proves the homogenization theorem. □ 


Remark 3.2. We remark that the same proof works in the case when (a t j) is not symmetric; 
indeed, it suffices to replace \ k above by the solution of the adjoint corrector equation. The same 
idea of proof can also be carried out in the case when (a^) are stationary ergodic random fields; 
indeed, the corrector equation in that case is much more involved but, by now, its solution and 
analogs of (13.41) and (13.51) . are well known. 


3.1. Decomposition of the homogenization error. To separate the fluctuations in the homog¬ 
enization error u £ — u that are due to the periodic oscillations in the diffusion coefficients from those 
due to the random potential, we introduce the function v £ which solves the following deterministic 
problem: 


“ SE (“«(?) = /M- 


x £ D, 


(3.7) 


v £ (x) = 0, x £ dD. 

Here, the potential field is already homogenized, and we expect that v £ — u filters out the effect of 
the random potential. The problem above is well posed and its solution v £ is given by Q £ f. 

The standard periodic homogenization theory yields that v £ converges weakly in i7 1 (D) and 
strongly in L 2 (D) to u, for any / £ H^ 1 (D). Using this function, we can write the homogenization 
error for ([1.11) as 

u £ — u = (u £ — v £ ) + ( v £ — u). (3.8) 

The deterministic part of the homogenization error is 

E u £ — u = E(u £ — v £ ) + ( v £ — u), (3.9) 

and the random fluctuation part of the homogenization error is 


u £ — Eu e = ( u £ — v £ ) — E(u e — v e ). 
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(3.10) 






The deterministic part of the homogenization error hence contains two parts, the mean of u £ — v £ 
and the periodic homogenization error v £ —u. Estimates for the second part amounts to convergence 
rate of periodic homogenization, and we recall some of the well known results below, together with 
uniform in e estimates of the Green’s function associated to Q s . We postpone the estimates for 
E(u £ — v £ ) to the next section. 

Theorem 3.3 (Estimates in periodic homogenization). Let D C be an open bounded C 1,1 
domain, v £ and u be the solutions to m and m respectively. Let G e (x,y), x,y € D, be 
the Green’s function associated to the Dirichlet problem of cud. Assume (A). Then there exists 
positive constant C, depending only on the universal parameters, such that 


(i) For any f <E L 2 (D), ||u £ - u \\ L 2 < Ce\\f\\ L 2 . 

(ii) For d > 2 and for any x,y G D, x / y, G e (x, y) satisfies 


C\x-y\ 2 ~ d , ifd+ 2 

C (1 + | log \x — 2 /||), ifd = 2 


and 

|VG e (*,y)| <C\x-y\ l ~ d . 


(3.11) 

(3.12) 


The 0(e) error estimates in L 2 was proved in [25] for d = 2, and in [15] for general O 1,1 domains; 
see also [20], The uniform in e estimates on the Green’s function and its gradient can be find in 
mm- The homogeneities in these bounds are the same as those for the Green’s function associated 
to constant coefficient equations, namely the Laplace equations. The striking fact that these bounds 
still hold for oscillatory equations is due to the fact that the problem (13.711 homogenizes. Periodicity 
or other structural assumptions on the coefficients are necessary. 


4. Estimates for the Homogenization Error 


In this section, we estimate the size of the homogenization error u £ — u. In view of the de¬ 
composition (13.811 , (13.91) , (13.101) and the error estimates in Theorem 13.31 it suffices to focus on the 
intermediate homogenization error u £ — v £ , with v £ = Q £ f defined in (|3.7|) . 

We introduce the function w £ which solves 


C £ w £ = -v £ v £ , 


(4.1) 


with homogeneous Dirichlet boundary condition. With the notations Q e and Q £,UJ introduced earlier, 
w £ is given by —Q £ v E v e . It follows that 

C £,u (u £ -v £ - w £ ) = -v e w £ , 

and u £ — v £ — u vanishes at the boundary. Hence we have u £ — v £ — w £ = —Q £,u i' £ w £ . Due to the 
assumption (A), Q £,u is uniformly (in e and u) bounded as a linear operator L 2 —> L 2 ; we have 

\\ u£ — vE \\l 2 < C\\w £ \\l2- (4.2) 


An estimate of u £ — v £ thus follows from the following result. 


Lemma 4.1. Let v £ = Q £ f and w £ as above. Under the same conditions of Theorem \2.A there 
exists a universal constant C and 


EIKII b(D) < 


Ce dM \\f\\ 2 L2 , 
Ce A \ loge|||/||| 2 , 


d¥= 4 , 

d = 4. 


(4.3) 


ii 




Proof. Using the Green’s function G e , we write 

w £ {x, UJ ) = J D G e{x, y)v v £ {y)dy. (4.4) 

The L 2 energy of w £ is then 

IKMIll, = J Ge(x,y)Ge{x,z)i> v ^ v £ (y)v £ (z)dydzdx. 

Take the expectation and use the definition of the auto-correlation function R of q, we have 

E||iV(-,u ;)||£ 2 = J D3 G £ {x,y)G £ (x,z)R v £ (y)v £ (z)dydzdx. (4.5) 

Integrate over x first. Apply the uniform estimates (13.111) and the fact (see e.g. Lemma A.l of 
my for any y 7 ^ z, 0 < a, (3 < d, 


dx 


ID 


| re — y\ d ~ a \x — z\ d ~P 


< 


We get 


ID 


\G e (x,y)G e (x,z)\dx < 


C, 

C(l + |log|y-«||), 
C\y - z\ a+ P~ d , 

'C\y - z\-^ d -^ A0 \ 
C(l + log|y- z\), 


Hence, if d > 2 and d 4, 

E|K(., W )|| \ 2 <c 


L 


\v £ (y)v £ (z)\ 


R 


y - z 


if a + f3 > d, 
if a + /3 = d, 
if a + f3 < d. 

if d± 4, 
if d = 4. 


dydzdx. 


(4.6) 


(4.7) 


£>2 | y — z\( d 4 ) v0 

When d = 4, the term (|y — z |( d - 4 )vo^-i s hould be replaced by 1 + | log \y — z\\. In any case, the 
above yields a bound of the form 


E\\w £ (-,uj)\\ 12 < C [ \v £ (y)\ (K £ * v £ ) (y)dy. 
JR d 


(4.8) 


Here, v e = v £ t d and Id denotes the indicator function of the set D, K £ (y ) = 7?(|)|y|^ 4 ~^ A0 ls p (y) 
if d / 4 and K £ (y) = i?(|)( 1 + lB p (y)\ log |y||) if d = 4. Here, B p is the ball centered at zero with 
radius p and p is the diameter of D. We check that, when d 7 ^ 4, 

f \R(y)\ 


\K‘ 


(y)\\p < [ 


R[- 


1 


;dy = 


f 

J 


= Ce 


dA4 


(4.9) 


|y|(^ 4)VO £-(d 4)VO |^|(^ 4)VO 

where in the last inequality we used R G L°° D L 1 (M d ). Similarly, when d = 4, 

\\K £ (y)\\ L i = [ ^r(-) (l + |log|y||)dy = e 4 / |-R(y)|(l + |log|ey||) < C7e: 4 |log£r|. (4.10) 

J B p Ve/ JB p/e 

To get the last inequality, we evaluate the integral on B\ and B p / e \ B\, and bound | log \ey\\ by 
| log \ep\\ for the second part. Apply Holder’s and then Young’s inequalities to (14.81) . we get 

E||^||| 2 <C||iL £ || Ll ||V||| 2 <G||iL £ || Ll ||/||| 2 . 

Combining this with the estimates in (14.91) and ()4. 101) . we complete the proof of the lemma. □ 
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4.1. Scaling of the energy in the random fluctuation. Now we estimate the L 2 (D) norm (the 
energy) of the random fluctuation u £ — E u £ which, in view of (13.101) , is the same as the fluctuation 
u e — v e — E(w £ — v £ ). 

Using the first order corrector w £ defined by (|4. 1 1) , and following the approach of H E3 , we 
represent u e — v £ in terms of a truncated Neumann series: 

u £ - v £ = -G e v e v £ + g £ v E g £ u £ v £ + g £ u £ g £ -u £ {u £ - v £ ). ( 4 . 11 ) 

In particular, the fluctuations in u £ — v £ can be written as 

u £ - e u £ = ~g £ v £ v £ + ( g £ v £ g £ v £ v £ - e g £ n £ g £ n £ v £ ) + ( g £ v £ g £ n £ (u £ - v £ ) - e g £ v £ g £ v £ (u £ - v £ )). 

The first term above is exactly w £ , which has mean zero and its energy was estimated in Lemma 
ECU The next lemma provides an estimate for energy of the second term in the above expansion. 


Lemma 4.2. Suppose that the assumptions of Theorem \2.h 1 are satisfied. Then there exists a 
constant C > 0, depending only on the universal parameters and /, such that 


ie II g £ v £ g £ v e v £ 


^g £ n £ g £ u E v £ \\ 2 L 2 {D) 


< 


Cs 2d 

Ce 8 1 loge| 2 
Ce 8 


if d = 2,3, 
if d = 4, 
if 5 < d < 7. 


(4.12) 


Let I| denote the left hand side of (14. 12|) : it has the expression 


if = E 


Id 5 


G £ (x, y)G £ {y , z) [v £ (y)n £ (z) - E n £ (y)n £ (z)\ v £ (z)dzdy ) dx 

Id \Jd 2 

G £ (x,y)G £ (x,y')G £ (y, z)G £ (y', z')v £ (z)v £ (z') 


y 


Eli/l-jvl — \ v[-)v\ — \ )-R 


y - z 


R 


y' - z' 


dz' dy' dzdydx. 


It is then evident that we need to estimate certain fourth order moments of ifix,ov), namely, the 
function 

^ u(x , y, t, s ) := E i fix)u(y)n(t)n(s) - (E v(x)v{y )){E u{t)u{s)). (4.13) 

Were v a Gaussian random field, its fourth order moments would decompose as a sum of products 
of pairs of R. This nice property does not hold for general random fields; however, the following 
estimate for p-mixing random fields provides almost the same convenience. 


Lemma 4.3. Suppose n(x,oj) is a random field with maximal correlation function g defined as in 
(12.61) . Then 

\y u (x,t,y,s)\ < d(\x - t\)d(\y - s|) + tf(|a; - s\)d{\y - f|) (4.14) 

where d(r) = (Kg(r/3))^, K = 4|M| L oo ( q X jD) . 


We refer to [18] for the proof of this lemma. Estimates of this type based on mixing property 
already appeared in [4]. We refer to [8] for an alternative way to control terms like and to 
section [3 for some comments on the connection of condition (S) with the lemma above. 
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Proof of Lemma \f.?\ Integrate over x in the expression of Ifj and apply the estimates (13.111) , (14.71) 
and (14.141) . We find, for d> 3, 


L 2 


(1 + ld= 4 | log | y - y'\\)\v £ (z)v £ (z')\ i y 


J d* \y — y'\^ d 4 ) v0 |y - z\ d 2 \y'— z'\ d 2 
(1 + ld= 4 | log | y - y'\\)\v £ (z)v £ (z')\ 


d 


+ 


d 


y - z 
£ 


d 


z — z 


z-y 


dz'dy'dzdy 


dz'dy'dzdy 


Id* \y - y'\ [d 4 l v0 |y — z\ d 2 \y'— z'\ d 2 

For d = 2, the terms | y — z | _ i rf_2 ) and | y' — z' | _ i rf_2 ) above should be replaced by 1 + | log \y — z\\ 
and 1 + | log \y' — z'\\ respectively. Let Ifi an< l I 22 denote the two terms on the right hand side of 
the estimate above. In the following, we set p to be the diameter of D. 


Estimate of Jfi • We use the change of variables 


y-y 


•-+ y, 


z — z 


i->- z, y' — z' 1 —^ y', z' i->- z'. 


£ £ 

Then the integral in becomes, for d > 3, 

^ r iyi z [ iy' f 


r(<f-4)V0 


' B 2 , 

p/e 


>B b 


Id |y| (d 4 l v 0 |y' + e(y - z)\ d 2 \y'\ d 2 


We integrate over y' first and apply (14.61) . then integrate over z' and obtain 

re /^ii„. £ || 2 _2d-2(d-4)V0 f (1 + ld=4| log M|)(l + lrf= 4 | log \e(y - z)\)d{y)d(z) ^ ^ 
\\v>e I „ \y\{d-i)vo\y _ z \(d-i)vo dydz - 


,B P/S 


When d = 3, the integral above is bounded because d E L 1 (M' i ) thanks to assumption (S), and 
we have Jfi + C£ 2d . When d = 2 , the situation is similar; after the integral over y' , there is again 
no singularity in the denominator. Hence, I £ n < C£ 2d . 

When d > 5, by Hardy-Littlewood-Sobolev inequality {251 Theorem 4.3], we have for p, r E 
(l,oo), 

f (w/i^yw 

Jr 2 d I y-z\ d 4 

Take p = and r = for any (d — 8 ) V 0 < 5 < d — 4. Then because d E L°° D L 1 (M d ), 
\y\ 4 ~ d E L P (B i), the above is finite and we have /|j < Ce 8 . 

When d = 4, we need to control the integral 

/ (1 + 1 log |ey||)(l + | log |e(y - z)||)d(y)d(z)dyd 2 ;, 

Jb 2 , 

p/e 

where D* = {y — ?/ — z + 2 7 | y, y', z, z' E D} is some bounded region formed by certain combinations 
of points in D. As a result, the logarithmic terms are bounded away from the poles. Hence, the 
above integral is bounded by 0 (| loge| 2 ), and If i + C + 8 1 logs] 2 . 




|y| 


d—4 


IL r 


Lp 


1 d-4 1 „ 

- H--h - — 2. 

p d r 


Estimate of I | 2 ■ We apply the change of variables 


y - z' 


*-+ y, 


y -z 


v. 


z' z'. 


£ 


£ 


14 


Z — z' !->• Z, 




























Then the integral in J 22 becomes, for d = 2, 


Ce 


2d 


' B 2 . 
p/e 


dydy' f dz' [ dz\v £ (z')v £ (z' + z)|(l + | log \z - ey||)(l + log \z - ey'\)d(y)d(y'). 
Jd Jb 0 


Integrate over z\ z and then over y' and y. We find that I 22 < Ce 2d . For d > 3, the same change 
of variables transforms J 22 to 

/ (1 + ld= 4 | log \z - e(y - y')\\)\v £ (z')v £ (z' + z)\ 


Ce 


2d 


/ dydy' dz 
J b 2 , J Bp Jd 

p/e P 


dz' 


d{y)d(y'), 


Id \z — e(y — y')\^ d 4 ) v0 |z — ey\ d 2 \z + ey'\ d 2 
After an integration over z' , we only need to control 

Ce™ f dydy' f dz- - + 7 K 

Jb 2 p/e Jb p \z-E\y-y )\ { 4,vu p — ey\ d 2 \z + ey'\ d 2 

When d = 3, an integration over z removes the singularities in the denominator. Then integrate 
over y and y l yields that / 22 < Ce 2d . 

When d > 5, we need to control the integral, after another change of variables e~ 1 z— (y — y 1 ) i— z 
and —y i->- y, 


dydy' dz- 


Ce 8 iJ(y)iJ{y') 


= dz 


\z\ d 4 \z — y\ d 2 \z — y’\ d 2 
with K(z) = {\y\~ ( ' d ~ 2 ' 1 * d)(z). Since d € L 1 n L°°(R d ), we have 


Ce 8 \K(z)f 
I ~\d—4 


\K(z)\ = 


d(y) dy 
d—2 


' bi ( z ) \z-y I 


+ 


/ 

JSL d \B 


d(y) dy 


< 


||$||l°° dy 


+ 


/ d(y)dy < C. 

jR d \Bx(z) 


d-4 d +4 

2d 


dz)\ z ~y\ d 2 Jb^z) \y - z\ d 2 

Moreover, by Hardy-Littlewood-Sobolev inequality, we have that 

\\K\\L 2 m = |||y|- (d - 2) *d(y)\\ L 2 {Rd) < C\\d\\ L2d/{d+ , )m < C\\d\\^\\d\\Ji 

Now we showed that K G L°° n L 2 (R d ). It follows that the integral to be controlled is finite and 
we have I 22 < Ce 8 . 

When d = 4, after the same change of variable as in the case of d > 5, we are left to control 


dydy' 


dz 


(1 + log \ez\ | )d{y)d{y')dydy'dz _ 8 


p/e 


/ B 


3 p/e 


| z — y\ 2 \z — y'\ 2 


= e 


7 (1 

J B 3p/e 


+ \log\ez\\)(K(z)) 2 dz, 


where K(z) = (lB p/e (y)\y\ 2 * d)(z). We verify again that K G L°° n L 2 for any 5 G (0,1). 

Estimate the integral again by breaking it into pieces inside and outside B\\ we find / 22 < Ce 8 \ loge|. 


Combining these estimates above, we proved (14.121) . 


□ 


Moving on to the last term in the series (14.111) . we observe that it cannot be controlled in the 
same manner as above. Indeed, the term u e — v £ is random and depends on i/(x,uj) in a nonlinear 
way. As a result, when we move the expectation into the integral representation, like in step (14.51) . 
we cannot get a simple closed form in terms of R. 

We hence choose not to address the interaction between u £ — v £ and the random fluctuation v £ 
in the potential directly. Instead, by an application of Minkowski’s inequality, we have 

||IE G e v £ Q £ v £ {u £ - u)\\ L 2( D ) < E|| g £ v £ Q £ v £ (u £ - u) || L 2 (D ). 
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Thus, we use the trivial bound on the L 4 (gi, L 2 (D)) norm of the fluctuations in g e v e g e i^{u £ — u): 

r e 2 := E\\g e i ' £ Q £ i> £ (u £ - u) - (Eg £ n £ g £ n £ (u £ - u)) || L 2 (D ) < 2E \\g £ v £ g £ v £ {u £ - u) || L 2 (D ), 

and only control the energy of the last term in (14.111) itself, as contrast to its variance. 

Lemma 4.4. Suppose that the assumptions of Theorem 12.31 are satisfied. Then there exists some 
constant C, depending only on the universal parameters and f, such that 


E\\g £ v £ g £ v £ {u £ -v £ )\\ L 2 {d] < 


Ce d 


if d = 2,3, 


Ce 4 1 logs1 2 if d = 4, 

if d > 5. 


(4.15) 




To prove this result, we estimate the operator norm of g e v e g £ , which is random since v £ depends 
on oj, and combine it with the control of u £ — v £ which is obtained already earlier. 

Lemma 4.5 (Mean value of the operator norm -*!?)• Under the same assumptions of 

Theorem \2.S\. there exists some universal constant C, such that 


e|| g £ u £ g 


e\\ L 2^ L 2 ^ 


< 


Ce d , 

Ce 4 \ loge| 2 , 
Ce 8 ~ d , 


if d = 2,3 
if d = 4, 
if d> 5. 


(4.16) 


Proof. For any h E L 2 (D), we have 


\\g £ n £ g £ h \\ 2 L2 = J yj ^ G E (x,y)u E (y)G E (y,z)h(z)dzdyj dx. 


Note that for almost every fixed iED, 

G £ (x, y)v £ (y)G e {y, z)h{z)dzdy 


>d 2 


It then follows that 


< W\\L 2 


[ G £ (x,y)u £ (y)G £ (y,-)dy 
JD 


L 2 


\\GeV e Ge\\ L *^i?(u) < J^yJ^G £ (x,y)is e (y,uj)G £ (y,z)dy) dzdx. 


Taking expectation, and we fined 
xn. n2 ^ f 

Jl 


mGeV e Ge\\L 2^ L 2 < 1^ G £ (x,y)G £ (x,r])R ^ ) G e (y, z)G E {y, z)dydydzdx. 


Integrate over z and x variables first. Using (13.111) and (14.61) . we find that the integrals over x and 
z variables are estimated as in (14.71) . Then we have 


e|| g £ v £ g 


en 112 < C 


e\\ L 2^ L 2 ^ 


l + Tc?=41 log \y y\\ 

r>2 \ \y — 77I 4 ) v0 


R 


y~v 


dydy. 


Change variables in the above integral and carry out the analysis as before. We find that (14.16[1 
holds. Note that the estimates become useless for d > 8. □ 
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Proof of Lemma \4-4\ For each uj E we have 

II g £ u e g £ u e {u e - v e )\\ L 2 < M\\G e v £ G e \\ L 2 ^ L 2 \\u £ - v £ \\ L 2 , 

where M is the uniform bound on the random potential in (12.41) . Take expectation and then the 
desired estimate follows from (14.161) . (14.21) and (14.31) . □ 

4.2. Scaling factor of the random fluctuations in weak topology. In this section we aim 
to find the correct scaling factor such that the random fluctuation u £ — E u £ , normalized properly 
according to this factor, converges with respect to the weak topology. For that purpose, we fix an 
arbitrary p> € L 2 (D) with unit norm, and estimate E ( u £ — E u £ , p) 2 . 

Using the series expansion formula (|4.11l) . we have 

(u £ - E u £ ,p) = -(G £ v £ v £ ,<p) + (l G £ v £ G £ v e v £ -E(G £ v £ G £ v £ v £ ),p) 

+ (G £ v £ G £ v £ (u E - v £ ) - E {G £ v £ G £ v £ {u £ - u £ )), p). 

Since the operators G e and G £ v £ G £ are self-adjoint on L 2 (D), we can move it to p. Set if e = G £ p. 
The above expression becomes 

(u £ -E u £ ,<p) = -{y £ v £ ,if £ ) + [(v £ G £ v £ v £ , if £ ) -E {y £ G £ v £ v £ ,if £ )\ 

+ \{v £ (u £ - v £ ),G e o £ ^ £ ) - E {v £ {u £ - v £ ),G £ v E ip e )\ (4.17) 
:=if + [Jf-EIf] + [Jf-EIf]. 

The aim now is to control the variances of If, j = 1, 2,3. 


Estimate for If. For If, which is mean zero, we have 
E(/[) 2 = e( / f v £ (x)v £ (x)ip £ (x)dx\ = f R 


J 

J R c 


>D 


Id 2 


x-y 

£ 


v £ (x)v £ ( y)if £ (x)if £ (y)dxdy 


[R £ * (v £ (y)ip £ (y)l D (y))\ (x)v £ (x)if £ (x)l D (x)dx 
^ C\\R £ \\L 1 {R d )\\ v£ 'lp £ ^D\\ 2 L 2(j SL dy 

Here, R £ (y ) = -R(|) is a shorthand notation. To obtain the last inequality, we applied Holder’s and 
Young’s inequalities. Note that |/? £ ||^i( R d) = e d ||i?|| L i( R d). Note also that /,</?€ L 2 (D) implies 
that v £ ,if £ € H 2 (D ) which is embedded in L 4 (D) for all 2 < d < 7. As a result, we conclude that 
E|Jf| < Cel 

Estimate for I\- Before calculating the variance of I s , we first confirm that E {If) 2 is larger than 
e d for d > 4. Hence, the target scaling £2 is correct (with respect to the weak topology) for the 
fluctuation u £ — v £ — E(ii £ — v £ ), not the difference itself. By direct computation, for d > 3, 


E(/f) 2 = 


R 


< 


!d a 


Id 4 


R 


x-y 

£ 

x-y 


R 


R 


x' — y' 


x' — y' 


G £ {x,y)G £ {x' ,y')v £ {y)v £ (y')^ £ {x)if £ {x')dx' dy' dxdy 


\v £ (y)v £ (y')if £ [x)if £ (x')| 
\x — y\ d ~ 2 \x' — y'\ d ~ 2 


(4.18) 


dy'dx'dydx. 


For d = 2, the last integral above should be replaced by 


' D 4 


R 


x-y 


R 


x — y 


v £ (y)v £ (y')ip £ ( x)iIj £ ( x')[ 1 + | log \x - y||][l + | log \x' - y'\ 


dy'dx'dydx. 
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After the change of variable 

x — y x' — y' . , . 

-h- x, - h-x, y^-y, y^-y, 

£ £ 

the integral to be controlled, for d > 3, becomes 


[ [ \R(x)R(x')\ 

Jb 2 , Jd 2 

p/e 


\v £ (y)i/> £ (y + £x)v e (y')ip £ (y 1 + £x')\ 

\ x \d-2\ x l\d-2 


dy'dydx'dx. 


Integrating over y and y' and then over x and x', we find that the integral above is finite. Hence, 
E(If) 2 is of order e 4 when d > 3. When d = 2, the change of variables in the logarithmic functions 
yield the term | loge| 2 , and we have E(If) 2 is of order e 4 | loge| 2 . We confirm that the scaling of the 
second term in (14. lip , when integrated against test functions, can be larger than or comparable to 
that of the leading term, when d > 4. 


We show next that if the mean is removed, the second term in (14.111) is smaller than the leading 
one in all dimensions. Using the definition of 4 a in (14.131) and the estimate in Lemma 14.31 we 
bound E (If — E/f) 2 = Var (If), for d > 3, by 


x y x' y' 


I <~TfUx, .yjc, ;</ yu yy yu yy yiy \xyip £lr ~ ''' , ™ A - ' 

£ £ £ £ 

'\ / y~y'\ |u e (y)u e (y')V ;£ (a;)V ,£ (a'’ / )l , /j / 


Var (If) = / G e (x,y)G e (x ,y )v £ (y)v £ (y )'ip £ (x)'ip £ (x)dx dy dxdy 

J D A \£ £ £ £ 1 

+ C f 


I® — y\ d 2 \ x ' ~ y'\ d 2 


-dx dy dxdy 


Id 4 


£ 


y - x'\ \^(y)^(i/W(x)r(x!)\ ix , d]Jixiy 


\ x _ y\d 2 |/J./ _ y!\d 2 


The second integral above is essentially the same with the first one if we interchange x' and y'. 
Hence, we focus only on the first one. After the change of variables 

x — x' y ~ y' / / / / 

- X, -^ y, y-x^-x, y^-y, 

£ £ 

the first integral becomes 


Ce 2d f dxdy f dx' f dy'd{x)d{y) 
lK 2d JB 0 JD 


| v e (y' + £y)v e (y')^ (y' — x' + £x + £y)il> £ (y' — x' + £y)\ 
\x' — £x\ d ~ 2 \x' — £y\ d ~ 2 


Integrate over y' hrst and use the fact that ||u e || i 4 < C and ||V ;£ ||l 4 (d) ^ C 1 - Then the above 
integral is bounded by, for d / 4, 



d(x)d(y) 


C£ 2d dx' 

x' — £x\ d ~ 2 \x' — £y\ d ~ 2 



Ge 2d ( d 4 ^ v0 d(x)i9(y) dx'dxdy 
\x — y|( d_4 ) v0 


where we integrated over x' and used (14.61) to have the inequality. The resulting integral is clearly 
finite. Hence we conclude that Var (If) < C£ 2d for 1 = 3 and that it is of order e d+4 for d > 5. 

When 1 = 2, there is only logarithmic singularity to start with in the expression of Var (If), and 
we find Var (If) < Ce 2d . 

When 1 = 4, the integral that remains after we integrate over x' has a term of the form (log |e(x — 
y)\)^-e(x-y)eB 2p - ^ follows then that Var (If) < Ce 8 ) loge|. 

To summarize, for d > 2, we have E|If — EIf| <C E|If|. That is, when the series expansion is 
integrated against test functions and the mean is removed, the second terms is much smaller than 
the leading term. 
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Estimate for J|. For the last term, we control it by the crude estimate E(/| — E/|) 2 < 2E(/3) 2 . 
From the expression If = (y £ (u £ — v £ ), Q E u £ if £ ), we have 

E |if | < E(\\is £ \\ l °o\\u £ - v e \\ L 2 \\Q £ u e ^ £ \\ L 2 ) < C (E\\u £ - v £ \\ 2 L 2 E\\Q £ v £ i/j £ \\ 2 l2 ) 5 . 

Since Q £ v £ ^ £ is exactly of the form of w £ defined in (|4.1[> . Owing to (|4.2I) and Lemma 14.11 we 
conclude that E |I|| is of order e d for d = 2,3, of order e 4 | log e| for d = 4, and of order e 4 for d > 5. 
Hence, for all 2 < d < 7, the truncation term in the Neumann series, with respect to the weak 
topology, has a scaling factor that is smaller than that of the leading term (which is of order £ 2 ). 

Remark 4.6. We find that for 2 < d < 7, the random fluctuation u £ — E u £ scales like £2 when 
integrated against test functions, and the leading term is the dominating one. We do not expect 
the dimension constraint d < 7 is not intrinsic. Firstly, it is related to the fact that we stopped at 
the second order iteration in the Neumann series, and had to control the last term by the crude 
estimate given by the Minkowski inequality (not taking advantage of removing the mean). Secondly, 
it is also needed when we claim that = Q £ v £ ip is in L 4 (D). In general, if we assume stronger 
condition, namely / £ C(D), then v £ is always bounded, and we only need if £ € L 2 (D) which holds 
in all dimensions if ip £ L 2 (D). 

We conclude this section by collecting the facts obtained above to give a proof of Theorem 12.31 

Proof of Theorem \2.3[ Let v £ be as defined in (13.711 . In view of (13.101) and Minkowski inequality, 
we have 

E||u £ - Eu £ || 2 2 < E [2||w £ - v £ \\ 2 L 2 + 2 ||E(u £ - u £ )||| 2 ] < 4E||u £ - u £ || 2 2 . 

Owing to (14.21) and Lemma 14.11 we have (12.161) . 

In view of (13.81) . (|4.2j) . Lemma [441 and Theorem 13.31 we have 

E||u £ — u\\ L 2 < E||u £ — v £ \\ L 2 + ||u £ — u\\ L 2 < CeH/11^2. 

This proves (12.151) . 

Finally, to estimate E|(u £ —Eu £ , (p)l 2 \ f° r an arbitrarily fied if £ L 2 (D), without loss of generality 
we can assume ||<^||l 2 = 1. Then this term is precisely what has been studied immediately above. 
With I £ , j = 1, 2, 3, defined earlier, we have showed that for 2 < d < 7, E | Yl^j=i(Ij ~ EZJ)| < Ce 2 , 
which is precisely (12.171) . □ 


5. Limiting Distribution of the Random Fluctuation 


In this section, we study the limiting distribution of the scaled random fluctuation, £“2 fu £ — E u £ ), 
in functional spaces. As mentioned earlier, the choice of space depends on dimension. When d = 1, 
convergence in law in C{D) of the random fluctuation was proved in |10l 4j. We prove Theorem 
IP below, which establishes convergence in law of the random fluctuation in L 2 (D) for d = 2,3, 
and in H^ 1 (D) for d = 4, 5. 

Multiplying £“2 to the series expansion (14.111) . we obtain the following expression for the scaled 
random fluctuation: 



(g £ u £ Q £ u £ v £ - e g £ v £ g £ v £ v £ ) (■ g £ v £ g £ v £ {u £ - v e ) - e g £ v £ g £ v £ (u e 



(5.1) 
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Our strategy, as in mi, is to prove that the leading term X £ = — e 2 Q £ v £ v £ contributes and 
converges in law to the right distribution depicted by Theorem 12.41 and show that the other terms 
converge in stronger mode to the zero function and hence has no contribution to the limiting law. 

At the purely formal level, all these steps are the same as in the setting of non-oscillatory 
diffusion coefficients. Indeed, we already established controls for the second and last terms above 
in the previous section. Moreover, the e dependence in Q e and v £ is not a problem, as we will see 
later, for the convergence of the characteristic functions of X £ , thanks to the fact that G e ip —>• Qip in 
L 2 for any tp E H~ 1 (D). This dependence, however, does impose difficulty on showing the tightness 
of the measures of {X £ } £ . As discussed in section [71 the old approach for tightness in [5] fails and 
new ideas are needed. 

Our new approach is to use some non-optimal but convenient tightness criteria, described in The¬ 
orems [AT] and [A21 for probability measures on H k (D) that are induced by processes in H k+S (D), 
k = —1,0 and s > 0. Since we do need s to be fractional in (0,1), we recall some definitions 
regarding fractional Sobolev space; see PU for reference. Given an open set K C M. d , the fractional 
Sobolev space Hq(K), for s E (0,1), is the closure of C^fi(K) in the norm 

11 112 _ 11 112 , f \ u ( x ) -u(y) I 2 , , 

\\ U \\h°(K) ■- \\ U \\l2(K) + J R2 | x _ y \d+2s dxd V- 
When K = an equivalent norm for u E H s (R d ) is 

/ a + ieHw 2 m- (5.2) 

v ' JR d 

Moreover, for s E (0,1), H~ S (K) is defined to be the dual space (Hq(K))', and in particular when 
K = the above norm for H~ s (R. d ) is still valid. 


5.1. Limiting distribution in L 2 (D) for dimension two and three. For d = 2,3, we prove 
that the leading term X £ in ()4.11l) converges in law in L 2 (D ) and show that the other terms vanish 
in the limit. The next lemma, together with Theorem lA.il yields tightness of X £ , which is the key 
step. 


Lemma 5.1. Suppose that the conditions of Theorem ] 2. .91 are satisfied. Assume further that d = 2, 3. 
Then for any s E (0, there exists a constant C, depending only on the universal parameters and 
s, such that 


E 


e 2g e i/ £ v £ 


2 

< C. 


H s 


(5.3) 


Proof. For each fixed uj E 11 and e > 0, e 2 Q e i> £ v £ belongs to Hq(D) and hence also Hq(D) for any 
s E (0,1). In particular, its H s semi-norm has the expression 


£ 2 G e V £ V £ 


J H S (D) 


ID 2 


\(G e v £ v £ ){x) - (g £ is £ v £ )(y)\< 

\rp _ ^y|c£-|-2S 


■dydx. 


Take expectation and use L 4 bounds of v £ , we have 


E 


'2Q e v £ v e 


2 

H S (D) 



1 (G £ (x,z) - G e (y,z))(G e (x,£) - G e (y,Q)\ 

\x — y\ d+2s 
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d^dzdydx. 
























Figure 1 . Decomposition criteria of the domain of integration based on the relative 
position between four points. Left: (x,y) E D 2 . Middle: (x,y) E DRight: 
{x,y) € D\. 





We claim: there exists G depending only on the universal parameters and s, such that for all 

f l(Q.(«.»)-g.(».»))(^.fl-0.fe.f))l dydx<c< (5.4) 

Jd 2 \x-y\ d+zs 

where the integral is understood in the principal value sense, as the limit of the integrals over 
D \ {B$(z) UBj(()}, 5 —>• 0. Since we obtain uniform in 5 bounds for these integrals, we simplify 
notation and denote the integration domain still by D 2 . We decompose the integration region D 2 
into three parts D 2 , j = 1,2,3 as follows: in D 2 , one of the points in {z,£}, namely £ without 
loss of generality, lies outside B p {x) U B p {y ) where p = \x — y\\ in D 2 , one of the points, namely z 
without loss of generality, lies in B p {x) and satisfies \z — x\ < \z — y\ and at the same time £ E B p (y) 
and |£ — y\ < |£ — x|; in D 2 , £ and z cluster around one of the points in {x,y}] without loss of 
generality, assume this point is x , so z,£ E B p (x) n {p \ \p — x\ < \p — y |}. In Figured] the relative 
positions between {x,y,z,f} are illustrated for each case. 

Let Ij be the integral over D 2 of the integrand in (15.41) . We estimate Ij, j = 1,2,3, separately 
and we focus on the case of d = 3. It is clear that when d = 2, the only change is that the Green’s 
function has logarithmic bound, and the analysis below can be adapted. 

On without loss of generality, we assume that \z — x\ < \z — y\ (if otherwise, we would switch 
the role of x and y). Hence | G e (x, z) — G e (y, z)\ < C\x — z\ 2 ~ d . By mean value theorem, 


|G e (a;,£) — G e (y,£)| < |VG £ (j 7 ,£)||x — y |, for some p between x and y. 

By the Gradient bound (13.121) and the fact that 1 77 — £| > \y — £|/2, we have 

\VG e (ji, £)| < ^ ^ hence \G e {x,£)-G e (y,£)\ < ^^d-i 

As a result, we have, 


h < 


G 


l D 2 \x — y\ d + 2s 1 \x — 


Id—2 


Is/ — ^1 


d—l 


dxdy. 


Integrate over x first and then over y. using (|4.6I) in each step; we hnd as long asO<s<^,/i<G 
for some G that only depend on the universal parameters and s. 
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On D%, we have \G e (x, z) - G e (y,z)\ < C\x - z\ 2 d and \G s (x, £) - G £ (y, £)| < C\y - £| 2 d . At 
the same time, \x — z\ < \x — y\ and \y — £| < \x — y |, so we may split the singularity into the 
integrals over x and y so that each of them is essentially not singular. That is, 

C 1 1 


I2 < 


L 


D% \x — z\ 2+S| 


\d-2 


dxdy. 


, al *-* r 3 |y-ei t+s \ x ~ z \ d ~ 2 \y~t\ c 

We note that the integral above can be separated, and as long as 0 < s < 2 — f = each integral 
is finite and hence I 2 < C. 

On D%, we assume without loss of generality that z and £ cluster around x. Then we have 
| G e (x,rj) — G e (y,rj)\ < C\x — y\ 2 ~ d for rj E {z,£}. At the same time, \x — y\ > \y — z |. As a result, 
we have 


h < 


L 


c 


1 


1 


dxdy. 


I D 2 I y — z\ d ~ r \x — z\ 2s+t \x — z\ d ~ 2 \x — C\ d ~ 2 
We choose r > 0 so the integral over y is uniformly bounded. The integral over x is also bounded 
as long as 2s + r < (4 — d) A 2 = 1, and we have I 3 < C. We note that for any s € (0, |) there 
exists r E (0,1 — 2s) satisfying the constraint 2s + r < 1. 

The above bounds are uniform in 5. Therefore, taking the limit 6 —> 0, we proved (15.41) . Integrate 
over z and £ in the integral expression of E [£~^Q £ i2 £ v e ] 2 Hs \ in particular, integrating R(-/e) yields 
a factor of e d that cancels the one in the denominator. We conclude that E [£~^G e v e v e ] 2 HS < C for 
each fixed s E (0, |). Combining this with E \\£~^Q e v e v £ \\ 2 L2 < C, which is due to (14.11) for d = 2, 3, 
we proved (15.40 . □ 


Remark 5.2. The key step in the proof above is to derive (15.41) . which concerns only the Green’s 
function G e and hence is obtained from a purely deterministic argument. Indeed, the scaling factor 

_d 

e 2 plays a role only afterward when we integrate against R £ , and it disappears in the final estimate 
because it is the right scaling for integrals of R e . In section [ 6 ] where we consider the case of long 
range correlated random potential q(x,u), the scaling in X s will be different, but the tightness of 
(the measures of) X £ , with the right scaling, is obtained in the same way as above. 


Next we address the convergence of the characteristic function of the measure of X £ . In view of 
Theorem 12.11 this amounts to proving 

Lemma 5.3. Assume (S). For any fixed tp E L 2 (D), we have 

— — ^=(G £ v £ v £ , <p ) L 2 d _ lstnbutl ° . n > j\f ( 0; a ^j ; where a 2 := a 2 J u 2 (x)(Gp>) 2 (x) dx. (5.5) 

Proof. Moving the operator Q e to (p, we have 

- i=(GeV £ v £ ,(fi) L 2 = - ]= [ v (-) v £ {x)'fi £ (x)dx 

Vs d Ve d Jd v £' 

where fi £ = G e tp. Let /f [tp\ denote the random variable above. Set if = Qtp and introduce 

Jf [ip\ := - -= [ v (—] u{x)if(x)dx. (5.6) 

V£ d JD v e y 

Since v(x,u) is a stationary ergodic random field that has short range correlation, mf E L 2 (D), we 
apply the well known functional central limit theorem, see e.g. Theorem in Bal [3], and obtain 

J{W\ — tribUtl ° - P > hW\ := cr (^J^G(x,y)u(y)dW(y),p S j ^ ~ AT ^0, a 2 j ^{u(y)if(y)) 2 dy^ . (5.7) 
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The last relation ~ above means equal in law. We note that 

E| v Q (v £ ^ £ - uif)dx^ < C\\v E r - utp f L 2 , 

and from periodic homogenization theory, we have v £ —>• u in L 2 , if £ —>• if in L 2 , as e —> 0; moreover, 
v £ and if £ are bounded in L°° since H 2 (D ) is embedded in L°°(D) for d = 2,3. As a consequence, 
the right hand side above converges to zero as e —>• 0. As a result, 

n M = J £ M + (IM - JIM) 

is a sum of a term that converges in distribution to I\[p\ and a term that converges to zero in 
L 2 (Pl). The desired result follows immediately. □ 

Finally, we collect the facts obtained above to give a proof of Theorem 12.41 (i). 

Proof of Theorem \2.4\ (i). Owing to Lemma 14.21 and Lemma 14.41 for d = 2,3, we have 

E \\e~i (g £ v £ g £ v £ v £ - Eg £ v £ g £ v £ v £ ) + e~i ( g £ v £ g £ v £ (u £ - v £ ) - E g £ v £ g £ v £ {u £ - v £ )) || L 2 < el 

By Chebyshev inequality, these two terms, as random elements of L 2 (D), converge in probability 
to the zero function. It follows that the limiting distribution of e - 2 (u £ — Em £ ) is given by that of 
the leading term X £ (u) := —e _ 2 g e u £ v £ . 

Let X be the right hand side of (12.181) . It is a random element of L 2 (D) defined on some 
probability space (O,^ 7 , P) on which the Wiener process W(y,u>) is defined. Let the distribution 
of X be P x and its characteristic function be cj) P ~ . We note that, for any tp € L 2 (D), the inner 
product (X,tp) has Gaussian distribution A^(0, cr 2 ), with <r 2 defined in (15.51) . Indeed, 

E pX (X, ip) = a E pX J (^j G(x, y)<p(x ) dx^j u(y) dW(y) = 0, 

and 

E pA (X, p) 2 = a 2 E pA ^ G(x, y)p(x)dx^j u{y) dW(y)^ = a 2 J^(gp) 2 u 2 dy, 

This shows (A, tp) ~ A/"(0, cr 2 ) in law. By Lemma f5.3l for any fixed p € L 2 (D), the random variable 
( X £ ,p ) converges in distribution to A/"(0,cr 2 ). This shows that, as mentioned in Remark 12.21 the 
characteristic function of the law of X £ converges to that of X. In view of Lemma [5. II and Theorem 
IA.1I the distribution of {A e } £g ( 0) i) in L 2 (D) is tight as well. Consequently, by applying Theorem 
EH we complete the proof of Theorem 12.41 (i). □ 

5.2. Limiting distribution in H~ l (D) for dimension four and five. For dimension d > 4, we 
do not expect e~^(u £ — E u £ ) to convergence in distribution in L 2 (D), because as shown in (12.161) . 
the fluctuations scale like f 2 | loge |2 for d = 4, and scale like e 2 for d> 5. In both cases, the scaling 
is much stronger than £ 2 . Nevertheless, we prove that convergence in law in H~ l (D) holds. 

As before, the key step is to show that the probability measure in H~ l (D) of the scaled leading 
term {X £ } := —e _ 2 g e i/ £ v £ in the expansion (14.111) is tight, and to show that the characteristic 
function of this measure converges. 
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Let us first address the characteristic function cj) X£ . We note that L 2 (D ) is naturally embedded 
to H~ 1 (D). For any / E L 2 , the linear form Lf : Hq(D) —>• R given by Lj( ip) = is clearly an 

element of iL _1 (D), and 


II-Mh-md) = su p L fM<\\f\\L2- 

M h i<1 

We henceforth identify Lf E H~ 1 (D) with / when / € L 2 (D). For any £ E H~ 1 (D), let l be the 
element in Hq(D) that is related to £ by Riesz representation. Then we have 

(f,£) H - HD) = L f (l) = (f,l). 

That is, the inner product of / E L 2 (D) with £ is the same as the L 2 inner product of / 

with l. As a result, Remark 12.21 applies for distribution on i/ -1 (D): to show (j) P converges to 
ip pX as characteristic function of distributions in it suffices to prove ( X £ ,h ) —>• (X,h) in 

distribution as random variables, for each fixed h E L 2 (D). 

Now we address the tightness of the measures of {X £ }. Our strategy is to control the mean of 
||X e ||#- S (£)) for some s E (0,1) and then apply Theorem IA.21 To this purpose, we first observe 
that X e E L 2 (D ) and hence X s E H~ S (D ) if we set 


: H°(D) 


by X £ (h ) = / X £ (x)h(x)dx. 


(5.8) 


ID 


For any h E Hq (D), the above clearly defines a continuous linear functional. Moreover, if we identify 
the function X £ with its extension to by zero outside D, the above also defines an element in 
H~ s (M d ). Since dD is regular (as a matter of fact, C 0,1 boundary is sufficient), any h E Hq(D) 
can be extended continuously to Eh E H s (W i ) which satisfies ||£7r||# S ( R d) < T\\h\\ H s^ D y, see [TT1 
Theorem 5.4]; by duality H~ s (M. d ) is continuously embedded in H~ S (D). If fact, we have 


I^IIh-»(D) : = sup (X £ ,w) l 2 = sup (X £ ,Ew) l 2 

™effg(D),|M|.H- S (.D)<i w£H s (d),\\w\\ H s( D) <i 

< sup {X £ ,v) L 2 < T\\X £ \\ H s 


(5.9) 


veH 


We note that ||^ e ||H- s (K d ) can be calculated using the formula (15.211 . 

Consider, for each fixed y E D, the Green’s function G e (-,y) for the Dirichlet problem (13.71) . 
Extend G e (-,y ) to R d by zero outside D, and let G\ denote the extended function. Then G\ defines 
naturally a linear form on H s (M. d ) by 


G\ 


H s ( R d ) 


h i— 


R, 

G v £ (h) := f G y e {x)h{x) dx = f G £ (x,y)h(x)dx, 

Jwi d Jd 


(5.10) 


provided the integral is finite. Since Q £ is self-adjoint and by the Green’s function representation, 
w(y) := G £ (h ) is the solution to the Dirichlet problem C £ w = h on D, with zero boundary condition. 
Note that the restriction of h on D is in H S (D). Invoking elliptic regularity, we find that w is 
bounded in H S+2 (D). Let s E (0,1) if d = 4 and s E ( 5 ,1) if d = 5, then by embedding theorem 
of fractional Sobolev spaces, H S+2 (D) C C°’ a (D) with a = s + 2 — | E (0,1); see JTBJ Theorem 
1.4.4.1], As a result, \G v £ (h)\ < Cj|/i||#s where C only depends on the universal constants and the 
index s. We hence proved the following fact: 
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Lemma 5.4. Assume (A) and q > 0. Identify G £ (-,y), for each fixed y € D, with the element in 
H~ s (R d ) defined above. Suppose s € (0,1) for d = 4 and s € (^, 1) for d = 5, then there exists 
C > 0, depending only on universal parameters and s, such that 

l|G , e('>y)ll,H'-' i (R<i) < C. (5-11) 


Using this fact and the Fourier transform formula for the H s (R rf ) norm, we can prove the 
following control of which, together with Theorem IA. 21 yields tightness of {X s }. 

Lemma 5.5. Suppose that the conditions of Theorem \2. A are satisfied. Assume further that d = 4, 5. 
Let sG (0,1) if d = 4 and s 6 (^,1) if d = 5. Then there exists a constant C > 0, depending only 
on the universal parameters and on s, such that 

E||X £ || 2 h -. (d) <C. (5.12) 


Proof. We identify X £ with the element in H S (D ) C H s (R d ) defined earlier. In view of (15.91) . 
we have 

11*11-00) < C||V'|| H - (Ri) = L / |^X«({)| 2 (1 + |{| 2 )-<if, 

e J R d 

where LT'X £ denotes the Fourier transform of the (extended) function X £ . Using the integral 
representation of X £ , we rewrite the above as 

\\ x£ \\h-s = “77 [ d£(l + |£| 2 ) _s / dxdy f dzdte^' ( ' x ~ y ^G £ (x,z)G £ (y,t)i' £ (z)v £ (z)i' £ (t)v e (t), 

£ JRd jR 2 rf .Id 2 

where the Green’s functions are extended by zero to M. d for their first variables. Take expectation 
in this formula; we have 


E \\X £ 


112 

ll/f-s 



y)G e (x,z)G e (y,t ) 

(l + l£l 2 ) s 


dxdydf 


) 


R 


z — t 


v £ (z)v £ (t)dtdz. 


We claim that for any z and t in D, 



e i ^-y)G £ (x,z)G £ (y,t) 

(i + iei 2 ) s 


dxdydf, 


< C. 


(5.13) 


Indeed, we recognize the quantity inside the absolute value sign to be 

/ <( I i^(oi 2 (i+i^- s 

JR C 


(1 + 1 ^ 


^(0 l 2 (i + iei 


2 \-s 


The term on the right hand side is precisely \\G z e || jy-s(Rd) HG* 11#-*^) • I n view of Lemma l5.4l we can 
apply (15.111) to get an upper bound for the quantity above and prove (15.111) . Then (15.131) follows, 
which in turn completes the proof. □ 


Finally, we conclude this section by collecting the facts above and prove Theorem 12.41 (ii). 

Proof of Theorem \2.4\ (ii). Step 1: Limiting distribution of the leading term. In view of Theorem 
IA.1I and Lemma 15.51 the probability measures on induced by {X £ } is tight. To check the 

limit of the characteristic functions of {P x }, it suffices to prove (12.141) . This is done in Lemma 
15.31 By Theorem 12.11 we conclude that X £ —>• X in distribution on H^ 1 (D) where X is defined to 
be the right hand side of (12.181) . 
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Step 2: Convergence to zero of the higher order terms. By Lemma 14.21 and d = 4, 5, we see 
that the second term in u £ — Ku £ , i.e. Q £ v £ Q £ v £ v £ — E Q e v £ Q e i/ e v £ converges, in L 2 (Ll,L 2 (D)) and 
hence in L 2 (Fl, H^ 1 (D)), to the zero function. Similarly, the remainder term Q £ v £ Q £ v £ {u £ — v £ ) — 
E Q £ n £ Q £ u £ (u £ — v £ ) converges to the zero function in L l (Q, H~ l (D)). These convergence results 
are stronger than the mode of convergence in distribution in H~ 1 (D). The proof of Theorem 12.41 
(ii) is thus complete. □ 


6. The Long Range Correlated Setting 


In this section, we consider the case when the random potential q{x,uf) is constructed from a 
long range correlated Gaussian random field, as described in the assumptions (L) in section [2j 


For the scaling of the homogenization error, we have the following analogue of Theorem 12.31 We 
focus on the error u £ — E u £ because, as seen earlier, the main contribution to the deterministic 
error E u £ — u comes from the periodic oscillation in the diffusion coefficients, and Theorem 13.31 (i) 
holds independent of the correlation length of v(x,uj). 


Theorem 6.1. Let D C be an open bounded C 1,1 domain, u £ and u be the solutions to (11.111 
and m respectively. Suppose that (A)(P) and (L) hold, f € L 2 (D). Then, there exists positive 
constant C, which depends only on the universal parameters, such that if 2 < d < 5 and 0 < a < d 
or 6 < d < 7 and 0 < a < 6, 


E ||u £ — Ku £ \\ L 2 < 


Moreover, for any tp G L 2 (D), 2 < d <7 and 0 < a < d, 



if 4, 
if d = 4. 


E \(u £ -Eu £ ,<p) L 2 \ < Ce 2 |M | L 2 ||/|| L 2 . 


( 6 . 1 ) 


( 6 . 2 ) 


The proof of (12.1611 is detailed in section 16.11 and the proof of (12.171) is found in section 16.21 This 
result shows that the random fluctuation u £ — Eti £ caused by the long range correlated random 
potential scales like e(" A4 )/ 2 in energy norm, and scales like e rf//2 with respect to the weak topology. 
Since a < d, we note that the random fluctuation in this setting is larger than the case of short 
range correlated potential. We mention that if a < 2, then the random fluctuations dominate the 
deterministic fluctuation caused by the periodic diffusion. 

The next result exhibits the limiting law of the rescaled random fluctuation e~^{u £ — Eu £ ). In 
the presentation, we define formally W a (dy) as W a (y)dy ; W a (y) is a centered stationary Gaussian 
random held with covariance function ~Ei{W a (x)W a {y )) = n\x — y\~ a , where E denotes the expec¬ 
tation with respect to the distribution of W a . Here, n = n g V 2 > 0, where n g and V\ are defined in 
(1271) and M . 

Theorem 6.2. Suppose that the assumptions in Theorem \6.1\ hold. Let k , be defined as in (USD 
and G(x,y) be the Green's function of (|1.2I) . Let W a (dy ) be defined as above. Then 

(1) For d = 2,3, as e —>• 0, 

U£ ~^ U£ } distribution > ^ [ G(x, y)u(y)W a (dy), in L 2 (D). (6.3) 

V£ q £ ->° Jd 

(2) For d = 4, 5, as e —>• 0, the above holds as convergence in law in H~ l {D). 
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Remark 6.3. The right hand side of (16.31) is an integral with respect to the multiparameter 
Gaussian random processes W a , we refer to [22] for the theory. Let X denote the result of the 
integral. When d < 4, a < 4, the Green’s function G(x, •) is in £ d /G-"/ 2 ) and X is a random element 
in L 2 (D). In general, X is understood through the Fourier transform of its distribution. Given 
h* € H~ 1 (D), 4> pX (h*) is defined to be Eexp(*-0« f D (G(-,y) 1 h*(-))u(y)W a (dy)). In particular, 
for any fixed positive integer N and functions {(pp, 1 < i < N} in L 2 (D ), the random variables 
Ii := {X,ipi) = y/R f D (G(-,y),ipi(-))u(y)W a (dy), i = 1, ■ ■ ■ ,N, is joint Gaussian, centered and 
have covariance matrix Ejj := E (Rlj) given by 


Eli .- K 


i 

Jd 2 


(uQ(pi)(y)(uQ(pj)(z) 


I V ~ 4 


dydz. 


(6.4) 


6.1. Scaling of the energy in the random fluctuation. In this section, we prove the first part 
of Theorem 16.11 which characterizes how does the random fluctuation m e -Em £ scale in the energy 
norm. 

The approach is the same as in the proof of Theorem 12.31 Lemma 16.41 below provides estimate 
for the leading term in the series expansion (14.IIP , and hence on E||u £ — v £ \\ 2 . This, together with 
the estimate in Lemma 16.51 on the operator norm, controls the remainder term of (14.IIP . Finally, 
Lemma 16.61 below controls the energy norm of the second term in (]4.1 1 [) . 

The first result below is already in [5]; we record the details here for the sake of completeness, 
and for demonstration of the key techniques. 


Lemma 6.4. Let Q e be as defined below (HU); let v(x,u) be the random, field constructed in (L). In 
particular, the autocorrelation function R(x) has the asymptotic behavior like n g \x\~ a with a < d. 
Then, there exists constant C that depends only on the universal parameters such that, for any 
f € L 2 (D), we have: 


E WGeVefWl* < < 


Ce aM \\f\\ 2 L2 , 
Ce a \\f\\l 2 , 

Ce 4 | loge|||/||^ 2 , 


if d 4 4, 

if d = 4, a < 4, 

if d = 4, a = 4. 


(6.5) 


Proof. The L 2 norm of G e n £ f has the following expression: 


\\GeV £ f \\ L 2 = J d yj D Gs(x, y)v £ (y)f(y)dy ) dx. 


After writing the integrand as a double integral and taking expectation, we have 

MGe^fWh = [ G e(x, y)G e (x, z)R £ (y - z)f(y)f(z)dydzdx. 

Jds 

Use (13.IIP to bound the Green’s functions. Integrate over x and apply (14.6P . We get 


nGs^ £ f\\h < c [ 

Jd 


1 + 1(2=41 log \y ~ z \ 


Id 2 \y - z\ G 4 ) v0 

Change variable (y, y — z) i->- (y, z). The above integral is bounded by 

f f 1 T 1(2=41 log \ z 

JdJb b k ~|( d - 4 ) v0 


I Re(y ~ z)f{y)f{z)\dydz. 


\R £ {z)f(y)f(y ~ 4\dydz. 


( 6 . 6 ) 


(6.7) 
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Here and in the sequel B p denotes the ball centered at zero with radius p which is the diameter of 
D. Integrate over y first, then we have: 


mGeV £ f\\ 2 L 2 < C\\f\\l, 


\Re(z)\{l + ld= 4 | log M 


-dz. 


( 6 . 8 ) 


Ib p M (d_4)v0 

Decompose the integration region into two parts: 

f D\ := < T} n B p , on which we have |i? e | < M 2 , 

\ D 2 := {| z£~ x \ > T} n B p , on which we have |i2 e | < Ce a \x\ 

where the bounds on R £ above are obtained by (L2) and by Lemma lA. 31 the constant T is defined 
in that lemma. The integration on D\ can be calculated explicitly, and it is bounded by 


(6.9) 


C 


S, 


1 T 1 d=41 log 121 


dz = C 


L 


eT r d 1 ^ _|_ | rf=4 | logr|) 

r (d-4)V0 


dr. 


IB eT 4)V ° 

When d 7 ^ 4, this integral is of order e dA4 . When d = 4, this integral is of order e 4 | loge|. 
Similarly, the integration over D 2 is bounded by 


Ce‘ 


I, 


1 + lrf=4| log |-g| 

,\B eT |^|( d - 4 ) v0 +“ 


dz = Cs ( 


L 


p 1 + lrf= 4 | log r\ 

1 g.rp T , (d —4)V0—d+l+Q! 


dr. 


When d 7 ^ 4, this integral is Cs a (p~ ( ' d ~ 4 ' >A0+d ~ a — an d hence is of order e aA4 . 

When d = 4, this integral has leading term Ce a (p 4 ~ a log p — (Te) 4_a | log \eT\ |), and hence is of 
order e a if a < 4 and of order e 4 | loge| if a = d = 4. 

Combining the above estimates on the integrals over D\ and D 2 , we obtain (16.5|h □ 


The above proof differs from that of Lemma l4.1l in the short range correlation setting only in the 
steps of controlling the integrals of R £ . In the long range correlation setting, because R(x) is not 
integrable, we cannot expect to gain a factor of e d by scaling the variable in R e . Instead, we gain 
a factor of e a by using the asymptotic of R £ outside an Te ball; see Lemma lA.31 This method of 
controlling the integrals of R £ will be used constantly in the sequel. 


Lemma 6.5. Under the same assumptions of Theorem \6.1l there exists some constant C depending 
only on the universal parameters, such that 




Ce a , 

if d = 2,3, 

Ce a , 

if d = 4, a < 4, 

Ce 4 | loge| 2 , 

Si. 

II 

S 

II 

C r e aA(8-d) > 

if d> 5. 


( 6 . 10 ) 


Proof. Following the proof of Lemma 14.51 we have 


n\Gev e G t 


ell L 2 ->L 2 — 


< c 


L 


1 + ld= 4 |log|y - 77 II 


r>2 \ | y — rj\( d 4 ) v0 

Proceeding as in the proof of Lemma 16.41 we obtain (16.101) . 


R 


y - v 


dydr\. 


□ 


In view of the relation (14.21) . we obtain immediately that, with u £ and v £ defined in (|l.lj) and 
(13.71) respectively, E ||u E — u £ || 2 2 satisfies exactly the same estimates in (|6.5I) . Combine this and the 
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above result, we find that ¥i\\Q e v £ Q s v £ (u e — v e )\\ L 2 is of order s a for cZ = 2,3,4 and a < 4, of order 
e 4 | log £ 1 2 if a = d = 4 , of order £ 2 +( aA ( 8 -< i ))/2 for d> 5. 

Lemma 6.6. Under the assumptions of Theorem \6.1l then there exists some constant C, depending 
only on the universal parameters, and 

E || Q e v e Q e v e v e - Eg e u £ g £ u £ v £ \\ 2 L 2 (D ) < £ a . (6.11) 


Proof. The expression for the quantity If := E ||C/ £ u £ t/ £ z7u £ — KQ £ n £ Q e i' £ v £ \\ 2 L 2 ^ D ' ) is found below 
(14.12[) . Integrate over x and apply the estimates (13. lip . (I4.7jl and (IA.10|) . We find, for d > 3, 7| 
can be bounded by 


E 

P#{(1,2),(3,4)}’ 


<7(1 + ld= 4 | log \y-y'\ 


d a I y — y'\^ d 4 ) v0 |y - z\ d 2 \y'— z 


— \R £ (x p(1) - x p(2 ))Re{x p ( 3) - a; p( 4))| dzdy'dzdy, 


t\d- 


where the ordered set {aq, X 2 , X 3 , X 4 } refers to {y, z, y',z'}, and p runs over the set U defined in 
(1A.9I) below, and p 7 ^ {(1,2), (3,4)}. When d = 2, the potential terms | y — z\ 2 ~ d and | y' — z '\ 2 ~ d , 
which result from an application of the bounds (|3 . 11 j) on the Green’s functions, should be replaced 
by | log | y — ,z|| and | log | y' — z'\\. In the summation above, there are 14 terms and each of them 
contains a product of two R £ functions evaluated at the difference vectors between a pair of points 
in {y, z, y ', z'}\ more importantly, at most one of the difference vectors is in {y — z,y' — z'}, in which 
case that R £ shares the same point of evaluation with one of the bounds on the Green’s functions. 


Figure 2. Difference vectors between four points. 



We can divide the fourteen choices of p into three categories as shown in Figure [2j In the 
first category, the two vectors in the correlation functions are linear independent with both of the 
vectors in the Green’s functions; in the second category, one of the Green’s function shares the 
same argument with one of the correlation function; finally in the third category, the vector in one 
of the Green’s function is a linear combination of the two vectors of the correlation functions. 


Case 1. In the first category, we consider a typical term of the form: 
C{ 1 + 1 < 2=4 1 log | y - y'\\)\v £ (z)v £ (z')\ 


21 Jd 4 |y — y'\^ d 4 ) v0 |y —z| d 2 \y'— z'\ d 2 

Change variable as follows: 


R £ (y — y')R e (z — z') \ dz'dy'dzdy. (6.12) 


y-y ^y, 


z-y' y', 
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2 — z' !->• z 


z' z'. 

















The integral can be bounded by 


,B3 P 


dy'dzdy / dz‘ 

Jd 


M 1 + ld=4| log \y\\)\v £ {z')v £ {z' + z)R £ (y)R £ (z)\ 


|^|(d 4)V0|y/ — y\d 2| y/ — £|d 2 

The integration over z' yields a bound by ||w £ || 2 2 . Then an integration over y' and an application 
of (14.61) yield 

(1 + td=i \log |y||)(l + 1 ^= 4 1 log | y - z\\)\R £ (y)R £ (z)\ 


I 21 ^ C\\V £ \\ 2 L2 f 
JBj 


|^|(rf 4)V0|^ — z |(^ 4)V0 


dzdy. 


(6.13) 


This integral is the same as in the estimate of Ifi i n the proof of Lemma 4.2. When d = 3, the above 
is essentially the square of the integral of R £ , and is of order e 2a if a < d and of order e 2d \ loge | 2 
if a = d. When d = 2, the integration over y' earlier eliminates the mild singularity posed by the 
logarithmic functions, and the same conclusion as d = 3 hold. 

When d > 5, applying the Hardy-Littlewood-Sobolev inequality again, we have for some 1 < 
p,r < oo 



\R £ (y)/\y\ d - 4 \\R £ (z)\ 


I y-z 


d -4 


dzdy < C 


R £ 

\y\d~4 


\\R £ \\l-{b p ), 

Lv{B p ) 


1 d -4 
_ 4- 

p d 


1 

H— — 2. 
r 


(6.14) 


If 0 < a < d, we can take p = d+4 j Q _ T , r = for some (4 — a) V 0 < r < d — a. These choices 
satisfy ar < d and p < Hence, according to Lemma [A.51 ||||are bounded by Ce a and 

\\R £ \y\~ d+A \W I s °f order o(l). If 5 < a = d < 7, we can choose p = and r = with some 
0 < r < 4. Then we check that \\R £ \y\~ d+A \\Lv ||-R £ ||z/- < Ce 8 <C e a . 

When d = 4, we note that the function K(x ) := l + ls 2p (x)| log |a;|| belongs to the weak Lebesgue 
space for all q € (l,oo). Apply the refined Young’s inequality [31 Theorem 1.5], we have 


[ I R e {y)K{y)K(y - z)R £ {z)\dzdy < C\\R e K\\ LP ( Bp) \\R £ \\ l -(b p ), ~ + ~ < 2 - (6-15) 

d B p p V 

When a < 4, we can choose p = r = with 0 < r < 4 — a. Then by (IA.6I) and (IA.5I) . 
||i2 £ A'||iP||i? e ||^ is of order e 2a . When a = d = 4, we can choose p = r = j— with 0 < r < 2. 
Then \\R £ K\\lp ||-R e ||L r is of order e 8 ~ 2r | loge| <C £ a . 


Case 2. In the second category, we consider a typical term of the form: 


T £ — 
1 22 ~ 


C(1 + ld= 4 | log \y ~ y'\\)\v £ (z)v £ (z')\ £ , jI,Uj 

in, /,,/ 1 (d— 4 )voL, ~irf— 21 ,,,/ y/id —2 \ R (y-y)R (y-z) I dzd v dzd v • ( 6 - 16 ) 


/d 4 |y — y'\^ d 4 l v0 |y — z\ d 2 \y'— z‘ 
This time we use the following change of variables, 


y ~y' i-ty, y' - ^ ^ y', y- z^> z , z' ^ z'. 


The integral can then be bounded by 

/C( 1 + 3Lrf=41 log 


,B3 P 


dy'dzdy / dz' 
Jd 


v £ (z')v £ (z' + y + y'~ z)\\R £ (y)R £ (z)\ 


\y\{d-A)V0\ z \d-2\yi\d-2 

Integrate over z' first and then integrate over y'. Then we get 

(1 + ld= 4 | log \y\\)\R £ (y)R £ (z)\ 


J-2.2 C 11 v £ 11 ^2 f 

Jsj 


|y|(d 4)VO|^|2 
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dzdy. 


(6.17) 





















When d = 3, the integral over y yields a term of order e a if a < 3, and of order e 3 | loge| if a = 3; 
the integral over z is of order e a if a < 2, of order e 2 | loge| if a = 2, and of order e 2 if 2 < a < 3. 
In all cases, /|2 ^ f° r d = 3. When d = 2, the singularity \z \~ d+2 is replaced by a logarithmic 
singularity, and the conclusion If 2 <C e“ holds. 

When 5 < d < 7, if a < 4, then the integral over y is of order e“; that over z is of order e 2A “ 
if a 7 ^ 2 and of order e 2 | loge| if a = 2. If a = 4, the integral over y is of order e 4 | loge| and that 
over z is of order e 2 . If a > 4, the integral over y is of order e 4 and that over 2 is of order e 2 . We 
conclude that if a < 6, then I |2 ^ £Q - 

When d = 4, the integral over y is of order e a if a < 4, and of order e 4 | loge| if a = 4. The 
integral over z is of order e a if a < 2, of order e 2 | loge| if a = 2, and of order e 2 if a = 2. In any 
case, I 22 ^ e“. 


Case 3. For the third category, we consider a typical term of the form: 
C{ 1 + lrf= 4 | log |y - y'\\)\v e (z)v e (z')\ 


T £ - 
J 23 ~ 


/d 4 |y — 4 l v0 |y~ z\ d 2 |y' — z'\ d 2 

Change variables according to 


R £ (y — y')R e {y — z')\ dz'dy'dzdy. (6.18) 


y-y H-y, 


y' y', 


y- z<-^ z, 


y — z' i->- z'. 


The integral can then be bounded by 

,C{ 1 + lrf= 4 | log \y\\)\v £ (y - z)v £ (y' + y - z')\\R £ (y)R £ (z')\ 


,B $ 


dy dz dz / dy 
J D 


|y|(d 4)V0|_r|rf 2| y - ~l \ d 2 

Integrate over y' and z and apply the ||u e ||i,oo bounds; then we are left to control 

(1 + 1^=41 log |y||)|i? £ (y)i?V)| 


Ih < Cll^lli 00 [ 

Jb* 


-dz dy. 


(6.19) 


|y|(^ 4 )V 0 |y — z'\ d 2 

This term can be estimated in a similar way as what we have done for (16.131) . The difference is that, 
in the application of Hardy-Littlewood-Sobolev inequality, the factor (d — 4)/d should be replaced 
by (d — 2)/d. In particular, for d > 3, we have, for p,r£ (1, oo), 


n 3 < c 


(l + lrf= 4 1 log \y\\)R £ 


|y|(d 4)VO 


L p (B P ) 


Ill'll L r (B p ), 


1 d -2 1 „ 

- H--t- 1 - - — 2 . 

p d r 


When d = 3, set p = r = |. Then the above is of order e 2 “ if a < |, of order e 5 | loge |3 if a = |, 
and of order e 5 if | < a < 3. In all cases, If 3 ^ £a - 

When 5 < d < 7 and 0 < a < 6 , we choose p = (/ - and r = —for some 0 < r < 6 — a, 

— — 7 1 a-\-Z—OL—T a+r 7 

and check that ra < d and p < This implies that ||i? £ ||ir- < Ce a , and ||i2 £ |y| _ci+4 ||LP is of 

order o(l). When a > 6 , we are only able to choose p, r satisfying the above constraint, such that 
ar > d and (d — 2 + a)p > d and I | 3 < Ce 6 . 

When d = 4, choose p = r = |. If a < 3, then ap = ar < 4 and both ||(1 + | log |y||)i? £ (y)||LP 
and 11 R £ 11 jjr are of order e Q , and the above quantity is bounded by Ce 2a . If a = 3, the above is of 
order e 2a \ loge| 2 . If 3 < a < 4, the above is of order £ 6 | loge|. In all cases, we have /| 3 <S £ a . 

When d = 2, the singular term in (I6.19P is K{y — z 1 ) with K{x) := ls 2 p (a;)| log |a:||. Note 
K € L2,(M d ) for all q € (l,oo), and hence /f 3 < ||I? £ ||i,p||i? £ ||L’' with any p, r € (l,oo) satisfying 
p~ 1 + r _1 < 2 . Set p = r > 1 so that p < 2 d/a, we find that ||i? e ||iP||i? e ||ir < Ce 2d / p <C e a . 
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Combining all of the results above, we finish the proof of the lemma. 

In view of the discussion at the beginning of this subsection, we proved (16.11) . 


□ 


6.2. Scaling factor of the random fluctuations in weak topology. In this section, we prove 
the second part of Theorem 16.11 which characterizes how does the random fluctuation u £ — E u £ 
scale when it is integrated against test functions. 

The approach is the same as in section f l .21 We fix an arbitrary tp € L 2 (D) with unit norm, 
and start with the expression (14. 17[) for ( u £ — E u £ ,ip). Let /}, j = 1,2,3, be defined as in (14.171) . 
Following the analysis in section I4~2l we get E (If) 2 < C||i7 £ ||^i(5 ), which is of order e a in view of 
a < d and (IA.5D . In the estimate of Var (If), (14.141) must be replaced by (IA.10D . so the analysis is 
more involved. We provide the details and allow a = d in the analysis of Var (If). 


The expression for the variance of /| is, for d > 3, 


Var (If) = / ) Ge{x,y)G e {x',y')v e {y)v £ {y')^ e {x)il) E {x , )dx'dy'dxdy 


' D 4 


< 


,?. L 


pH( 1 , 2 ),( 3 , 4 )} 


C\v £ {y)v £ {z')^ £ {x)il) £ {x')\^ f , , 

\ x - y\d- 2 \ x , - y,\d ~2 \ R e(x p (i) ~ x m )R e {x p{3) - x m )\d y dx dydx , 


where in the summation p runs over U \ {(1,2), (3,4)}, and the ordered set {xi,X 2 ,x 3 ,x^} refers 
to {x,y, x',y'}. For d = 2, the above bounds \x — y\~ d+2 and \x' — y'\~ d+2 should be replaced by 
1 + | log \x — y\\ and 1 + | log \x' — y'\\. These integrals in the sum can be estimated in the same 
way as those of in the proof of Lemma 16.61 by considering three categories of p. We note that the 
integrals here are less singular than before because the absence of |y — y'| term. 


Case 1. In the first category, we consider a typical term of the form: 


T £ — 

J 91 — 


L 


21 Ip 4 \x — y\ d 2 \x' — y 

Change variable as follows: 


_ x , w{y _ y)| dyUxUydx , 


( 6 . 20 ) 


x — x' i->- x, y — x' x', y — y' i->- y, y' i->- y'. 


The integral can be bounded by 

,C\v £ (y')v £ (y' + y)'t/j £ (y' + y — x')i/) £ (y' + y + x — x')R £ [x)R £ {y)\ 


/ dx'dydx / dy' 
J Bl Jd 


\x' — x\ d 2 \x' — y\ d 2 


The integration over y' yields a bound by C , ||u e ||^ 4 ||V’ e || 2 4 . Then an integration over x' and an 
application of (14.61) yield 

(1 + lrf= 4 | log \x - y\\)\R £ (x)R £ (y)\ 


ih < c 


f 

J Bp 


■dydx. 


( 6 . 21 ) 


| X _ y|(rf-4)V0 

When d = 3, the above is essentially the square of the integral of R £ , which is of order e 2a if 
a < d and of order e 2 “|loge| 2 if a = d. Hence, /fx <S e“. When d = 2, the integration over x' 
earlier eliminates the mild singularity of logarithmic functions, and the same conclusion holds. 

When d > 5, applying the Hardy-Littlewood-Sobolev inequality, we have for some 1 < p, r < oo 
\R £ (x)R £ {y)\ . . __ 1 d- 4 1 


! B 2 \x - y \ d ” 4 dydx - C ^ R£ ^ Lp ( b p) W R£ h r (Bp)i 


-h 

p 


d 


+ - = 2 . 

r 


( 6 . 22 ) 
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If 0 < a < d, we can choose p = d+4 _f Q _ T , r = for some (4 — a) V 0 < r < d — a. Then 
ar < d and hence ||.R e ||i/- < Ce a . On the other hand, ||I? e ||i> is of order o(l). If a = d. we 
choose p = -Ap and r = with 0 < r < d — 4. Then ar > d and ap > d. It follows that 
\\R £ \\lp\\R £ \W < Ce d(p 1+r ' 1 which is of order e d+i . In all cases, we have Ifi ^ £ a - 

When d = 4, the \x — y \~ d+4 above is replaced by K(x — y) with K := 1 + 1b 2 (x)| log |x||, and 
K belongs to L„,( W l ) for all q € (l,oo). By refined Young’s inequality, we have 

J I R £ {y)K(y - z)R £ {z) I dzdy < C\\R £ \\ L p ( b p )\\R £ \\l-(b p ), \ + \ <2 ' ( 6 - 23 ) 

If a < 4, we can choose p = r = with 0 < t < 4 — a. Then ra < d and ||i? e ||i,p||I? e ||i,r- is of 
order e 2a . When a = 4, we can choose p = r = : with 0 < r < 2. Then is of 

order s s ~ 2t . In all cases, I £ n <C e". 


Case 2. In the second category, we consider a typical term of the form: 
f C\v £ (y)v £ (y')i/j £ (x)iJj £ (x')\^ 

Id 4 


J 22 ~ 


I® — 2 |x' — y'\ d 2 


■|i? £ (x — x')R £ (x — y)| dy'dx'dydx. 


(6.24) 


This time we use the following change of variables, 

x-yrp-x, x — x'x', y i->- y, 


x' — y' t-p- y'. 


The integral can then be bounded by 

, f C\v £ (y)v £ (x + y — x' — y')ip £ (y + x)ip £ (x + y — x')R £ (x)R £ (x')\ 
dx dx / dy - ; — , , „ -. 

Jd \x\ d ~ 2 \y \ d ~ 2 

Integrate over y first and then integrate over y'. Then we get 

I e 22 < CW^WlooWWl, [ \E^)E^l dx'dx. (6.25) 

J B 2 \ x \ 

The integrals over x' and x are separated. When d > 3 and a < d, the integral over x' is of order 
e a , and that over x is always of order o(l). Hence /| 2 £ a ■ When a = d, the integral over x' is 
of order e d \ loge|, and that over x is of order at most e 2 . In all cases, / 22 <C £ a . 

When d = 2, the above quantity amounts to ||i? e | log |x|| \\l 1 (b p ) \\R £ \\l 1 (b p ) which is of order e 2a 
if a < 2 and of order e 4 | loge| 3 . In all cases, / 22 -C £ a ■ 



Case 3. For the third category, we consider a typical term of the form: 

r C\v £ (y)v £ (y , )'ip £ (x)'il> £ (x')\. .. , ... . . 

d 2 !/ rid— 2 \ R ^ X - x )R (x - y ) I dy'dx'dydx. 


T £ — 
1 23 ~ 


! d a \x - y\ d 2 \ x' - y 
Change variables according to 


x — x' x', x — y (->• y, x — y y . 


The integral can then be bounded by 

C\v e {x — y)v £ {x — y')^ (x)^ (x — x')R £ (x')R £ (y') \ 


Bp 


dy dx dy J dx- 


\y\d 21 ^./ _ yl\d 2 


Integrate over x first and then over y; we obtain 

^23 ^ ll^lli 2 


£ii 2 , 1-112 / \R £ (x')R £ (y')\ 


B 2 \x' — y'\ d 2 


-dx dy 


(6.26) 


(6.27) 
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Applying Hardy-Littlewood-Sobolev inequality again to the integral above, we have, for d > 3 and 
for p, r € (1, oo), 

J 23 < C \\ R£ \\lp(B p ) \\ Re \\L-{B p ), ~ H + - = 2 - 

When d > 3 and 0 < a < d, set p = d+2 d a _ T and r = with (2 — a) V 0 < t < d — a. 
Then ar < d and p > 1. Hence ||i? e ||x,r is of order e Q and ||i? e ||i> is of order o(l). If a = d, 
then set r = and p = with 0 < r < d — 2. Then ra > d and pa > d , and hence 

11 R £ 11 lp 11 R£ 11 L r < Ce d ^~ 1+r ~^ = Ce d+2 . In all case, Jf 3 < e Q . 

When d = 2, the singular term in (|6.27l) is K(x' — y') with K{x) = 1 B 2p ( x )\ log Ml- Note that 
K € for all q € (l,oo), and hence I 23 < ||-R £ ||lp|M £ ||l t ' for any p, r € (l,oo) satisfying 

p~ 1 + r -1 < 2. Set p = r > 1 so that p < 2d/a, we hnd that ||i? e ||LP||-R £ ||i,r < Ce 2d / p <C £ a . 

Combining all of the results above, we proved that there exists some constant C that depends 
only on the universal parameters and the data ip, such that for any d > 2, 0 < a < d, 

Var (If) < Ce a . 

Finally, for the truncation term If, we have E|/f| < C (E||rt £ — u e ||^ 2 5 which, in 

view of (16.51) . is of order e“ A4 or £ 4 | loge|, much smaller than s '2 for all d <7. 

We hence conclude with (16.21) . In particular, for 2 < d < 7, the leading term in the series 
expansion (14.111) is dominating in the weak topology, and it scales like £ 2 . 


6.3. Limiting distribution of m £ -Em £ . We multiply the series expansion (14.111) with the proper 
scaling factor of the long range correlation setting, and get the following expression for the scaled 
random fluctuation e” ( u £ — Ek £ ): 


G e v e v e , (G £ v £ G £ v £ v £ -EG £ v £ G £ v £ v £ ) (G £ h> £ G £ v £ (u £ - v £ )-¥*G £ v £ G £ v £ (u £ - v £ )) 


+ 


+ 


(6.28) 


Following the approach of section (5j when 2 < d < 5 and 0 < a < d, the first term above dominates 
and this term alone contributes to the limiting distribution in L 2 (D) for d = 2,3, and in H^ 1 (D) 

a 

for d = 4,5. We hence set X £ = —e _ 2 'Q £ b £ v £ . 


Tightness of the probability measures of {X £ }. In view of Remark 15.21 the tightness of the prob¬ 
ability measures induced by {W} on L 2 (D ), when d = 2,3, follows from Lemma 15.11 with e _ 2 
replaced by e”. Similarly, the tightness of the probability measures induced by {X £ } on H~ 1 (D), 
when d = 4,5, follows from Lemma 15.51 


Limits of the characteristic functions of {X £ }. We have shown in section [5] that for the limit of 
the characteristic functions of the measures of {X £ }, whether in L 2 (D) for d = 2,3 or in H~ 1 (D) 
for d = 4, 5, it suffices to characterize the limiting distribution, for an arbitrarily fixed ip € L 2 (D), 
of If \ip\ := (X £ ,p) = ~£~ 2 (i/ £ v £ ,if £ ), where f> £ = Q e <p. Following the proof of Lemma 15.31 
we introduce Jf[p\ = —£~^ (v £ u,if) with = Qp. For the random field u(x,uj) = <h(g(x, w)) 
constructed in (L), applying Lemma 4.3 in [5], we have 

Ji[<p] = — 4 - [ v (-) u(y)fj(y)dy dlstnbutlo . n > f u (y)f(y) W a (dy ), (6.29) 

£2 J D \£l e^O J D 
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as convergence in distribution of real random variables. Since If [<p\ — [<p\ converges in L 2 (fl) to 

zero, Ii[<p\ has the same limiting distribution. This establishes the convergence of characteristic 
functions of the measures of X £ . 

Proof of Theorem \6.Sl The arguments above proved that X s converges in law to X, in L 2 (D) for 
d = 2, 3 and in H~ l (D) for d = 4, 5, where X is defined as the right hand side of (16.311 . In view of 
(16.111) . the second term in (16.281) converges in L 2 (fl, L 2 (D)) to the zero function. Similarly, due to 
(14.21) . (16.51) and (16.101) . the last term in (16.281) converges in L 1 (ll, L 2 (D)) to the zero function; see 
the remark below (|6. 101) . The conclusions of Theorem 16.21 follows immediately. 

7. Further Discussions 

In this section, we make some comments on the assumptions of the random potentials used in 
this paper, and discuss the case in higher dimensions d > 6. 

An alternative condition for (S). In the short range correlation setting for v(x,uj), we 
assumed the condition (S). Upon applying Lemma 14.31 we can bound the (partial) fourth order 
moment ^> u by the sum of two terms, each consisting product of a pair of functions 1 ? € L°°nL 1 (W i ). 
However, as remarked earlier, (S) essentially requires the mixing coefficient g(r), and hence i?(|.x|), 
to behave like o(r~ 2d ) at infinity, much stronger than R(x) being integrable. 

We remark that (S) is assumed mainly to simplify the presentation and the o{r~ 2d ) decay of q 
is not necessary. In fact, an alternative assumption used in [S] to control fourth order moments 
is: there exists 1 ? : ^ M + in L 1 D L°°(M d ), such that (1 A. 101) hold. This is clearly a much more 

general assumption, and it is satisfied if v(x,oS) = ${g{x,u)) with <1 satisfying (L2) and (L3), and 
g(x,u>) being a centered stationary Gaussian random field with correlation function R g = o(|x| -d ) 
as |x| —>• 00 . 

The conclusions of Theorem 12.41 still hold if (S) is replaced by the above alternative assumption. 
Indeed, we only need to modify the control of in the proof of Lemma 14.21 and in section 14.21 
For instance, in the first inequality in the proof of Lemma 14.21 we have more but finitely many 
integrals instead of two on the right hand side. Nevertheless, in all of these integrals, at most one 
of the function d has the same variable with one of the Green’s function, and all of them can be 
controlled in the same spirit of the proof of Lemma 16.61 

Comparison with the case of non-oscillatory diffusion. The main results of this paper 
shows that the framework developed in m 00 , in the setting of non-oscillatory differential operator 
with oscillatory random potential, applies even when the differential operator is also oscillatory, 
as long as we have uniform in £ control of the Green’s functions and their gradients, i.e. (13.111) 
and (|3.12|) . and provided that there is no random correlation between the diffusion coefficients 
and the potential. At the formal level, there is no difference in the proof, and the usual strategy 
using (truncated) series expansion applies. However, the role played by the oscillatory diffusion 
coefficients becomes prominent in getting the tightness of the measures of {A £ }. 

Let us recall the previous method used for tightness in the setting of non-oscillatory differential 
operator. Set C := — Yltj=\ A; Qxdx + and consider the Dirichlet problem (£ + v £ )u £ = / in 
D with zero boundary condition. Then u £ homogenizes to u, the solution of (11.21) . As in [5], the 

limiting distribution of e~^(u £ — Eu £ ), say, in the short range correlation setting, is characterize 
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by that of X £ = —e ~2 Qv £ u. To prove tightness of the measures {X £ } in L 2 (D), the strategy of [5] 
is to use the spectral representation of L 2 (D). Note that £ is formally self-adjoint and its inverse, 
i.e. Q, is compact on L 2 (D). Hence, £ admits real eigenvalues {Afc}^ =1 such that, 

0 < q < Ai < A2 < • • • , Afc —>• 00 as k —> 00 , 

and eigenfunctions {4>k} < kLi, ||0/c||l 2 = such that 

j — A k (f> k , in D, 

1 4>k = 0, on dD. 

Moreover, {4> k } form an orthonormal basis of L 2 (D) and we have the following representation of the 
space T~L°(D) = L 2 (D) and the Sobolev space P l {D ) = Hq(D)-, see Section 6.5]: for s = 0,1, 


n a (D) = lfeC°°(D) : ^(/,^)i 2 A|<oo , and \\v\\ 2 w := Y (/, K- (7-1) 

l k =1 J k=l 

A natural criterion for tightness of (the measures of) {X £ } is that their measures do not scatter to 
higher and higher modes. More precisely, let Pj y denote the projection operator in L 2 (D) to the 
space Wj\T := span {<j> 1 , • • • , 4>n} spanned by the first N modes. Then {X £ } is tight if E||Af e || i 2 < C 
and 

lim sup E\\X £ - P N X £ \\ L 2 = 0. (7.2) 

N ~>°° £6(0,1) 

Using the representation formula in (|7.1h . and the fact that Q(p k = (X k )~ l (j) k , we have 

1 OO 1 OO .. 

nx £ -p N x £ f L2 = - d y n^ e u^ k ) 2 = Td y 

£ k=N+l £ k=N+l 

As in section [4~2l we have sup eg ( 0)1 ) sup fc E(i/ £ u, 4 > k ) 2 < C. In view of the Weyl’s asymptotic formula 

2 

for the eigenvalues, X k ~ k<* for k large, we conclude that 

OO OO 1 

sup E||V-P W V||| 2 < £ T 7 < £ TT- 

ee (°’ 1 ) k=N+l k=N+l k d 

Hence, for d = 2, 3, we obtain tightness of {X £ } for free, as byproduct of the analysis in 14.21 
In the setting of this paper, £ above is replaced by C £ defined in (13.11) . The above approach 
for tightness fails completely. On the one hand, if we replace the eigen pairs (A k ,4>k)k by (A|, (p k ) k 
where the latter solve the eigenvalue problems associated to £ e , then instead of (17.21) . we obtain 

lim sup E||X £ — P^X S \\ L 2 = 0, 

N^oo £e (o,l) 

where Pfj is the projection to W ^ := span{^»f, • • • , <jf N }. This is useless because, a priori, the basis 
)k changes with e, and it is not clear that the union (over e € (0,1)) of unit balls in Wfj is 
still compact for all N. On the other hand, if we fix a spectral representation, say, using (X k ,4>k)k 
defined before. Then we no longer have the relation Q £ 4> k = (A k )~ l (j) k . It is not difficult to check 
that ||Vt/e0fc|| i 2 « -T= and this estimate is sharp. An application of Poincare inequality yields 
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that 11 £/£</>& 11 £2 < C/sfXk ~ k d , with C uniform in e and k. It is not clear at all how to improve 
this estimate. Consequently, in view of the estimate on If in section 14.21 we have 

1 oo oo oo oo 

sup E||A £ — P n X £ \\ 2 L 2 = — V E(z ' e u,g £ (j > k ) 2 < V CWQMl^ V V —.. 

^(0,1) e fc= V +1 k ^ +1 k^ + l Xk k^+l^ 

This fails to show (17.21) or the tightness of {X £ }, even for d = 2. 

In view of the analysis above, we find that the above approach for tightness, which is natural 
for non-oscillatory differential operators, fails completely in the presence of fast oscillations in the 
diffusion coefficients. The new approach used in sections 15. II and 15.21 is necessary and more stable. 

Higher dimensional case. We comment on the restrictions of dimensions in our results. Let us 
consider the short range correlation setting only. With respect to the weak topology, we expect that 

_ d 

the scaling factor for the random fluctuation should be e 2 in all dimensions. More precisely, for 
any fixed ip E L 2 (D ) 1 we expect that e _2 (u £ — E u £ , ip) converges in distribution for all dimensions. 
However, in this paper we prove this only for d < 7; see the control of Var(/|) in section 14.21 
This constraint is used to control the truncation term in the series expansion (14.111) . and we do 
not expect it to be intrinsic. In fact, if we have more information about the random field, we can 
truncate later in (14.171) until the last term is under control, and use higher order moments estimate 
in the spirit of Lemma 14.31 and Lemma lA. 51 to control the terms in between. 

When certain norm topology is considered, however, we expect that the dimension plays an 

_d 

intrinsic role in the choice of the functional space. Indeed, for the leading term X £ = —e 2 Q e u £ v £ 
to converge in law in L 2 (D ), it is necessary that E||A £ ||^ 2 is controlled uniformly in e. Since the 
homogeneity \x — y\~ d+ 2 of the Green’s function is relatively more and more singular in higher 
dimensions, this is simply not true. In the setting of this paper, we expect that convergence in 
law in H~ k (D), for certain k > 0 increasing with respect to d could be proved. To establish this 
in detail, we need again to push the series expansion further and have more information of the 
random held. We note that other choices of functional space could be convenient for this endeavor 
as well, in particular, the space defined in (17.11) for larger s, or certain negative Holder spaces; 
we refer to [SI El for such considerations. 


Appendix A. Some technical results 

A.l. Tightness criteria for probability measures in functional spaces. We hist present a 
tightness criterion for the probability measures {P XE } eg ( 0) i) on L 2 (D) induced by {A £ (-,w)} that 
are random elements in HfffD) C L 2 (D), s € (0,1]. 

Theorem A.l (Tightness in L 2 (D)). Let {A e (-, w)} eg (o,i) be a family of random fields on the 
probability space (fl, J 7 ,P), X £ {-,oj) € Hq(D) for some 0 < s < 1, for each fixed e E (0,1) and 
oj € fL Suppose there exists C > 0, independent of £ and ui, such that 

E ||A £ ||h s < C. (A.l) 

Then the family of probability measures {-P A " e } £ e(o,i) on L 2 {D) is tight. 
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Proof. By assumption P A? concentrates on the subspace Hq(D). For any fixed 5 > 0, set Mg = 
C(5 _1 and define 

As = {f £ Hq(D) | \\f\\ H s < M s }. 

Clearly, As is closed and bounded in Hq(D). In light of the fact that the embedding H^(D) M- 
L 2 (D) is compact [28] . we note that As is a compact set of L 2 (D). Now for any fixed e E (0,1), 
applying Chebyshev inequality, we find 


P x \As) = n{X £ €HZ(D), 


> 1 - 


E||X £ 


\H S 


Mg 


> 1 - 


X £ \\ h * < M s }) 
C_ 

M~s 


= 1 - 6 . 


l-P({||X e ||^ > M S }) 


Since 6 and e are arbitrary, the above shows that {P A " E } £ g(o,i) is tight. 


□ 


Next we give a similar tightness criterion for probability measures {P x ’}ee(o,i) on H 1 (D) 
induced by {X e (-,w)} which belongs to a smoother space. 

Theorem A.2 (Tightness in P[~ l (D)). Let {X £ (-, w)} e g(o,i) be a family of random fields on the 
probability space (n,P,P), X e {-,u) E H~ S (D) for some 0 < s < 1, for each fixed e E (0,1) and 
u € fh Suppose there exists a constant C > 0, independent of e and c o, such that 

E ||X e ||^_ 3 < C. (A.2) 

Then the probability measures {P A "~} £ g( 0 , 1 ) on H^ 1 (D) is tight. 


Proof. Since D is a bounded open set with regular boundary, the embedding H^(D) ^ Hq(D), for 
any 0 < s < 1, is compact SSI Theorem 1.4.3.2], By duality, the embedding H S (D) H 1 {D) 
is also compact. The rest of the proof is exactly the same as in the proof of the theorem above. □ 


A.2. Functions of long range correlated Gaussian random field. Here we record some 
results for the random potential v(x,oj) = $(g(x,u)) that is constructed in (L). In particular, we 
express the asymptotic behavior of its correlation function R(x), and derive a (partial) fourth order 
moment for v. 


A.2.1. Autocorrelation function of the long range model. 

Lemma A.3. Assume (L1)(L2) and let v(x,oj) be as constructed there. Set V\ = E{<7o < h(5'o)}) 9x 
being the underlying Gaussian random field in (L). Then there exist constants T, C > 0, depending 
only on the universal parameters, such that the autocorrelation function R(x) of q satisfies 

| R(x) - V?Rg(x)\ < CR 2 (x), for all \x\ > T, (A.3) 

where R g is the correlation function of g. Further, 

\^{g(y)q(y + z)} - ViRg(x)\ < CR 2 g (x), for all |x| > T. (A.4) 

The proof of this result can be found in [6, [5]. It says that v{x,oS) inherits the heavy tail from 
the underlying Gaussian random field. The next result describes estimates on the integrals of R, 
possibly against some potential function that has singularity at the origin. 
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Lemma A.4. Assume (L1)(L2) and let u(x,oj) be as constructed there. Let R(x) be the correlation 
function of u, R £ (x) be the rescaled function R{e~ l x), and B p be the open ball centered at zero with 
radius p. Then there exists C > 0, depending only on the universal parameters, on p and p below, 
such that 


(i) For any 1 < p < oo, 

\\ r£ \\lp(b p ) < < 

(ii) For d = 2,4, and for any 1 < p < oo, 

\\R £ {y)\log \y\\\\ LP{Bp) < 


Ce a , 


ap < d, 


Ce a \\oge\p, ap = d, 


Cep, 


ap > d. 


m 


For d > 3, and for any 1 < p < 


Ce a , 


R £ 


\y\ 


d—2 




LP(B P ) 


(iv) 


For d > 5, and for any 1 < p < 


Ce a , 


R £ 


\y\ 


d- 4 


< 


LP(B P ) 


(A.5) 


Ce a , 

ap < d , 


Ce a | 

loge 1+ p, ap = d, 

(A.6) 

d 

Cep\ 

log e|, ap > d. 



ap < d — p(d — 2), 


1 

ge\ p > 

ap = d — p(d — 2), 

(A.7) 

1 -2) 

i 

ap > d — p(d — 2). 



ap < d — p(d — 4), 


l 

g^K 

ap = d — p(d — 4), 

(A.8) 

l- 4) 

ap > d — p(d — 4). 



Proof. The proof of (i) can be found in [5j. We present the proof of (iv) here; the other terms above 
can be proved in exactly the same way. Using the decomposition (|6.9I) . we control the integral of 
\Rt\P\y\-v{d-4) on d 1 anc [ jp 2 separately. 

On D i, we have 



\R £ \p 

|y|p( rf - 4 ) 


dy < || R £ 


up 

IIl°° 



1 

r p( d ~ 4) 


dr < C£ d - p{d ~ A) . 


On D 2 , we have 

r \R £ \ P d < £Qp f p r d ~ l ^ = f Ce ap ( r d-«P-p(d-4)) \ P T£ } ap^d- p(d - 4), 

lB p \B Te \y\ p (d " 4) V ~ £ Jts rP( rf - 4 )+“P r [C£ ap (log r)\^, £ , ap = d - p{d - 4). 


Combining these estimates, we proved (1A.8D . 
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A.2.2. Fourth-order moments of i/(x,ui). Finally, we present a non-asymptotic estimate for the 
four-moments of n(x,oj) constructed in (L1)(L2), with the additional assumption (L3). In the 
following, we denote by U the collections of two pairs of unordered numbers in the set {1,2,3,4}, 

U ■= {p = {(p(l),p(2)), (p(3),p(4))} | p(i) € {1,2,3,4}, p( 1) ^p(2), p(3) ^p(4)}. (A.9) 

As members in a set, the pairs (p(l),p(2)) and (p(3),p(4)) are required to be distinct; however, 
the two pairs can have one common index. There are three elements in U that collect all four 
numbers. They are precisely {(1, 2), (3,4)}, {(1, 3), (2,4)} and {(1,4), (2,3)}. Let U * denote the 
subset formed by these three elements, and let U* be its complement. 

Lemma A.5. Assume (L) and let v(x,uj) be as constructed there. Then there exists d : M. d —>• R+, 
bounded and satisfying d(x) ~ |a;| _ “ as |x| —>• oo, and some C > 0 depending only on the universal 
parameters, such that for any four points {xi € M d ; 1 < i < 4}, 

4 

«n^)- z R(x p ( 1 ) - X p { 2 ))R{x p $) -Xp( 4 )) < C ^ &(x p (i)-X p (2))'&(x p (3)-X p (4)). (A.10) 

i =1 p&U, p&U* 

We refer to [8J Proposition 4.1] for the proof of this result. In particular, t? above can be chosen 
as the autocorrelation function R(x) of u(x,ui). As discussed earlier, (lA.lOj) can be viewed as an 
alternative for the estimates in Lemma 14.31 
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